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r^, Abstract 

We perform an extensive study of FCNC and CP Violation within Supcisymmet- 
^ ric (SUSY) theories with particular emphasis put on processes governed by 6 — > s 

transitions and of their correlations with processes governed by 6 — )■ d transitions, 
^ s d transitions, D'^ — oscillations, lepton flavour violating decays, electric 

^ dipole moments and [g — 2)^. We first perform a comprehensive model-independent 

(3 analysis of AF = 2 observables and we emphasize the usefulness of the -Rfc — 7 plane 

^ in exhibiting transparently various tensions in the present UT analyses. Secondly, 

we consider a number of SUSY models: the general MSSM, a flavour blind MSSM, 
the MSSM with Minimal Flavour Violation as well as SUSY flavour models based on 
^ abelian and non-abelian flavour symmetries that show representative flavour struc- 

^ tures in the soft SUSY breaking terms. We show how the characteristic patterns of 

correlations among the considered flavour observables allow to distinguish between 
these different SUSY scenarios. Of particular importance are the correlations be- 
tween the CP asymmetry 8^4, and Bg — > between the anomalies in S(f,Ks ^-^d 
between S^Ks ^^'^ de, between S^^ and {g — 2)^ and also those involving lep- 
ton flavour violating decays. In our analysis, the presence of right-handed currents 
and of the double Higgs penguin contributions to Bg mixing plays a very important 
role. We propose a "DNA-Flavour Test" of NP models including Supersymmetry, 
the Littlest Higgs model with T-parity and the Randall-Sundrum model with cus- 
todial protection, with the aim of showing a tool to distinguish between these NP 
scenarios, once additional data on flavour changing processes become available. As a 
byproduct, we present the SM prediction for BR{B+ t+u) = (0.80 ± 0.12) x 10"^ 
that follows solely from an analytical formula for this branching ratio in terms of 
AMs,d and S^Ks asymmetry and which does not involve Vub and Fb uncertainties. 
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1 Introduction and Motivation 



The Standard Model (SM) of elementary particles has been very successfully tested at the 
loop level both in flavour-conserving electroweak (EW) physics at the LEP and the SLC and 
also in low-energy flavour physics. 

On the other hand, it is a common belief that the SM has to be regarded as an effective field 
theory, valid up to some still undetermined cut-off scale A above the EW scale. Theoretical 
arguments based on a natural solution of the hierarchy problem suggest that A should not 
exceed a few TeV. 

Besides the direct search for New Physics (NP) at the TeV scale (the so-called high- 
energy frontier) that will be performed at the upcoming LHC, a complementary and equally 
important tool to shed light on NP is provided by high-precision low-energy experiments (the 
so-called high-intensity frontier) . In particular, the latter allows to indirectly probe very short 
distances even beyond those accessible by direct detection. 

In the last years, the two B factories have established that the measured Bd flavour- 
and CP-violating processes are well described by the SM theory up to an accuracy of the 
(10 — 20)% level [H |2]. Unfortunately, irreducible hadronic uncertainties and the overall good 
agreement of flavour data with the SM predictions still prevent any conclusive evidence for 
NP effects in the quark sector. 

This immediately implies a tension between the solution of the hierarchy problem that 
requires NP at the TeV scale and the explanation of the Flavour Physics data in which this 
NP did not show up convincingly. An elegant way to avoid this tension is provided by the 
Minimal Flavour Violation (MFV) hypothesis O IH El El E], where flavour violation, even 
beyond the SM, is still entirely described by the CKM matrix. As a result, it turns out that 
a low energy NP scale at the level of few TeV is still fully compatible with the flavour data 
within this minimalistic scenario [8],[9]. 

In this context, the question we intend to address in this work is whether it is still possible 
(and to which extent) to expect NP phenomena to appear in the Bg system where the SM has 
not been experimentally tested at the same accuracy level as in the 5^ system. In particular, it 
is well known that 6 — )■ s transitions represent a special ground where to perform efficient tests 
of NP scenarios tl0llIIlll2l[l3l[I3lll5lll6llI7lll8]. Indeed, CP violation in 6 ^ s transitions is 
predicted to be very small in the SM, thus, any experimental evidence for sizable CP violating 
effects in the Bg system would unambiguously point towards a NP evidence. Recent messages 
from the CDF and DO experiments seem to indicate that this indeed could be the 

case [211122]. 

On the theoretical side, there exist many well motivated NP scenarios predicting large 
effects especially in 6 — )■ s transitions. Among them are supersymmetric (SUSY) flavour 
models based on abelian [23l[22[25l|26l|271[28l[29l|30lini[3^ non- 
abelian [3Hll31linillIlll21ll31ll3llSllinill3llHlllSEni[^ flavour symmetries naturally 
leading to large NP contributions in 6 — )■ s processes while maintaining, at the same time, 
the NP contributions occurring in s — )■ c? and (sometimes) b ^ d transitions, under control. 
Moreover, also Grand Unified Theories (GUTs) represent a suitable ground where large NP 
effects in 6 — )■ s transitions can be generated [TOl [TTl |53l ISH |55l [561 Ell- fact, GUTs link 
leptonic and hadronic sources of flavour and CP violation and the observed large atmospheric 
neutrino mixing is transmitted to a large flavour violation in 6 -H- s transitions fTDl fTTj . 

In the present work, we focus on the NP predictions for the CP violating and CP conserving 
6 — )■ s transitions within SUSY models and their correlations with other observables, measured 
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in coming years in dedicated flavour precision experiments in K, B, D and charged lepton 
decays. In particular, when we deal with specific models, the source of the fiavour and CP 
violation for 6 — )■ s transitions will simultaneously generate not only AB = 2 and AB = 1 
processes that will turn out to be correlated, but will also have impact on observables outside 
the B meson system. The major aim of the present work is twofold: 

i) to quantify the NP room left for 6 — )■ s transitions compatible with all the available 
experimental data on AF = 2 and AF = 1 processes, 

ii) to outline strategies to disentangle different NP scenarios by means of a correlated 
analysis of low energy observables, including also K and D systems as well as lepton 
fiavour violation, electric dipole moments (EDMs) and the anomalous magnetic moment 
of the muon {{g — 2)^). 

Since many analyses along this subject appeared in the literature [2], [121 lEl EHl [59], |60l 
l56t [61] , we want to emphasize here the novelties of our study. In particular, 

• We consider a very complete set of AF = 2, AF = 1 and AF = processes, i.e. the 
EDMs and the {g — 2)^. To best of our knowledge, the current analysis represents the 
most complete analysis present in the literature in this subject as for i) the number 
of processes considered and ii) the inclusion of all relevant SUSY contributions. We 
believe that the steps i) and ii) are extremely important in order to try to understand 
the pattern of deviations from the SM predictions that we could obtain by means of 
future experimental measurements. 

• We perform an updated analysis of the bounds on the fiavour violating terms in the 
SUSY soft sector, the so-called Mass Insertions (Mis), in the general MSSM, in the light 
of all the available experimental data in fiavour physics. 

• In addition to the general MSSM, we study several well motivated and predictive SUSY 
frameworks, such as the MSSM with MFV, a fiavour blind MSSM (FBMSSM) and in 
particular SUSY fiavour models based on abelian and non-abelian fiavour symmetries 
and compare their predictions with those found in the Littlest Higgs Model with T- 
parity (LHT) and a Randall-Sundrum (RS) model with custodial protection. 

• We outline a comprehensive set of strategies to disentangle among different NP scenarios 
by means of their predicted patterns of deviations from the SM expectations in many 
low energy processes. At the same time, we exploit the synergy and interplay of low 
energy flavour data with the high energy data, obtainable at the LHC, to unveil the 
kind of NP model that will emerge, if any. 

Our paper is organized as follows. In sec. |2j we present a model independent analysis of 
Bd, and B^ mixing observables, analyzing in particular the existing tensions between the 
data and the SM in the i?^ — 7 plane. We also investigate what room for NP is still present 
in the Bd and Bs systems. 

In sec. |3} we list the formulae for AF = 0, 1, 2 processes in a general MSSM that are most 
relevant for our purposes. 

In sec. |4| we briefly review basic features of SUSY models with abelian and non-abelian 
flavour symmetries as well as models respecting the MFV principle. We also address the 
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question of the stability of the squark mass matrix textures, predicted by the flavour models, 
under RGE effects from the GUT scale to the EW scale. 

As our paper involves many observables calculated in several supersymmetric models, 
we outline in sec. |5] our strategy for the numerical analysis that consists of five steps. The 
first step is a model independent analysis within the general MSSM framework resulting in 
the allowed ranges of the Mis. The next two steps deal with three specific supersymmetric 
flavour models in which right-handed currents play an important role: the abelian model by 
Agashe and Carone (AC) based on a U{1) flavour symmetry and the non-abelian models 
by Ross, Velasco-Sevilla and Vives (RVV) ^T] (or, more precisely, a specific example of the 
RVV model, i.e. the RVV2 model [62]) and by Antusch, King and Malinsky (AKM) [52] 
based on the flavour symmetry SU{3). These three models are then compared in the last 
two steps with a flavour model with left-handed (CKM-like) currents only [l2l |63], with a 
flavour blind MSSM (FBMSSM) [Ml ESI ESI EH EH] and with the MFV MSSM with additional 
CP phases [SSI EOl EI]- The numerical results of our strategy are presented systematically 
in sec. |6| The correlations between different observables, offering very powerful means to 
distinguish between various models, play an important role in this presentation. 

In sec. [7], we review very briefly the patterns of flavour violation in the Littlest Higgs Model 
with T-Parity (LHT) [72l[I71[73l[71l[Z5l[76l[77] and a Randall-Sundrum (RS) model with 
custodial protection [TJl US] and compare them with the ones identified in supersymmetric 
models in sec. |6| The main result of this section is a table summarizing and comparing the 
sensitivity of various observables to NP effects present in each model. This table can be 
considered as a "DNA-Flavour Test" of the considered extensions of the SM. 

In sec. [8| we exploit the complementarity and the synergy between flavour and LHC data 
in shedding light on the NP scenario that is at work. 

Finally, in sec. |9| the main results of our paper are summarized and an outlook for coming 
years is given. A compendium of one loop functions can be found in appendix |A] and in 
appendix |B] some details on the convention of the MSSM parameters used throughout our 
paper are shown. 



2 Model-independent Analysis of AF = 2 Observables 
2.1 The Bd System 

The present unitarity triangle (UT) analyses are dominated by AF = 2 processes. We begin 
by reviewing the status of the UT trying to outline transparently possible hints of NP and 
related tests to falsify or to confirm them. 

We consider the following two sets of fundamental parameters related to the CKM matrix 
and to the unitarity triangle: 
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are the length of the sides of the UT opposite to the angles j3 and 7, respectively. The latter 
are defined as follows 

Vtd = \Vu\e-'^ , Kb = iKftle-*^ , (2.5) 

and one has 

sin 2/3 = ^^^(| - ^ tan 7 = ^ . (2.6) 
+ (1 - qY q 



The parameter set in (2.1 ) can be obtained entirely from tree level decays, hence, it should 
be unaffected by any significant NP pollution (see however |78]). The corresponding UT is 
known as the reference unitarity triangle (RUT) ^79j. In contrast, Rt and /3 in the parameter 



set in (2.2) can only be extracted from loop-induced FCNC processes and hence are poten- 
tially sensitive to NP effects. Consequently, the corresponding UT, the universal unitarity 
triangle (UUT) [3j of models with constrained minimal flavour violation (CMFV) IHIHQ], could 
differ from the RUT signaling NP effects beyond not only the SM but also beyond CMFV 
models. Thus a comparative UT analysis performed by means of these two independent sets 
of parameters may represent a powerful tool to unveil NP effects. The dictionary between 
these two sets of variables is given by 



Rh = 



-—2 — — ^ 1- Rt cos l3 

' Rf - 2Rt cos /3 , cot 7 = „ ■ ^ , (2.7) 

Rt sm (j 



Rt = ^/l + i^2-2i^5Cos7 , cot /3 = ^ ^^^"'^ . (2.8) 
V Rf, sm 7 

Assuming no NP, the parameters Rt and (3 can be related directly to observables 



where AMd and AMg are the mass differences in the neutral Bd and Bg systems, S^^Xg 
represents the mixing-induced CP asymmetry in the decay Bd — t- ipKs and the value of the 
non-perturbative parameter ^ is given in tab. |3| In the presence of NP however, these relations 
are modified and one finds 



Rt = d,IW^J^J^ , sm{2P + 2<Ps,) = S^Ks, (2.10) 



A V ^B, V AMs V C 



where AMg = AM^^Cb^ and (pB^ is a NP phase in Bd mixing defined analogously to (2.18) 
below. 

The experimentally measured values of AMd and S^Ks can be found in tab. [l] along with 
their SM predictions. While in AMd, there is still room for a NP contribution at the 25% 
level, scenarios with large new CP violating phases are strongly constrained by the bound on 
SipKs- This is also illustrated by fig. [T| showing the experimental constraints in the complex 
Mf^ and the S^Ks-^^d/ AMf^ plane, respectively. 

2.2 The K System 

In the K system the central role at present is played by e^, that represents another crucial 
ingredient of any UT analysis. 
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observable 


experiment 


SM prediction 


exp./SM 


AMk 
ex 


(5.292 ± 0.009) X lO'^ ps"^ EI] 
(2.229 ±0.010) X 10-3 [81] 


(1.91 ±0.30) X 10-3 


1.17±0.18 


AM, 

SipKs 


(0.507 ±0.005) ps-i P 
0.672 ± 0.023 [IJ 


(0.51 ±0.13) ps-i 
0.734 ± 0.038 


0.99 ±0.25 
0.92 ±0.06 


AM, 
AMrf/AM, 


(17.77 ±0.12) ps-i (82| 
(2.85 ±0.03) X 10-2 


(18.3 ±5.1) ps-i 
(2.85 ±0.38) X 10-2 


0.97 ±0.27 
1.00 ±0.13 



Table 1: Experimental values and SM predictions for AF = 2 observables. The SM predictions 
are obtained using CKM parameters from the NP UTfit [83]. The last column shows the ratio 
of the measured value and the SM prediction, signaling the room left for NP effects in the corre- 
sponding observable. We do not give a SM prediction for AMk because of unknown long distance 
contributions. 



The CP-violating parameter ex can be written in the SM as follows 

e^|SM = n^C^BK\Vcb\'^\Vus\^ (]-\Vcb\'^R^sm2^r]uSo{xt) + Rt sin ^{r]ctSo{xc,Xt) - ii^cXc 



2 

where is a numerical constant 



(2.11) 



and the SM loop function 5*0 depends on Xj = mj{mi)/M^ (where rriiimi) is the quark 
mass rrii computed at the scale rrii in the MS scheme) and can be found in appendix |a| The 
factors ?7tt, 77^ and r^cc are QCD corrections known at the NLO level [HU [851 EHl [87], Bk is a 
non-perturbative parameter and is explained below. 

While eK has been dormant for some time due to a large error in the relevant nonper- 
turbative parameter and CKM parameter uncertainties, the improved value of B^-, the 
improved determinations of the elements of the CKM matrix and in particular the inclusion 
of additional corrections to ex [HB] that were neglected in the past enhanced the role of this 
CP-violating parameter in the search for NP. 

Indeed it has been recently stressed [88] that the SM prediction of ex implied by the 
measured value of sin 2/3 may be too small to agree with experiment. The main reasons for 
this are on the one hand a decreased value of Bk = 0.724 ± 0.008 ± 0.028 [HH] (see also [UP]), 
lower by 5-10% with respect to the values used in existing UT fits [831 EI]! and on the other 
hand the decreased value of in the SM arising from a multiplicative factor, estimated as 

= 0.92 ±0.02 [HHip. 

Given that eK oc Bxt^e, the total suppression of eK compared to the commonly used for- 



mulae is typically of order 15%. Using directly (2.9) together with (2.11), one finds then 

|e^|SM ^ Y8 ^ 0.25) x 10-^ , (2.13) 
to be compared with the experimental measurement [81] 

le^l^^p = (2.229 ± 0.010) x 10-^ . (2.14) 



^Using instead the values of g and 77 from tab. [sjone finds the SM prediction in tab. [ij 



6 



2 Re(Mf2)/AM^'^ 
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Figure 1: Left: The complex plane. The thick solid circles show the constraint on the 
mixing amplitude from the measurement of AM^/AM^ assuming no NP in AMg, while the thick 
dashed circles correspond to the constraint coming from AM^ alone. The thin lines correspond to 
NP contributions to with the indicated fixed phases 6^, defined analogously to (2.18). The 
white region corresponds to the allowed region of S^Ks (see tab. [l]) and the SM value corresponds 
to the black point. Right: Same as in the left panel, but now in the S'^i^-g-AM^/AM^'^ plane. 



/cb\'^ and R^. 



The 15% error in (2.13) arises from the three main sources of uncertainty that are still Bk, 

Wry 



On general grounds, the agreement between (2.13) and (2.14) improves for higher values 
of Bki Rt or \Vcb\ and also the correlation between ex and sin 2/3 within the SM is highly 
sensitive to these parameters. Consequently improved determinations of all these parameters 
is very desirable in order to find out whether NP is at work in S^Ks or in ex or both. Some 
ideas can be found in [93], EH [921 |9l] . We will now have a closer look at possible tensions in 
the UT analysis in a somewhat different manner than done in the existing literature. 



2.3 A New Look at Various Tensions: -R5-7 Plane 

While the above tension can be fully analyzed by means of the standard UT analysis, we 
find it more transparent to use the i?6 — 7 plane [HSl [M] for this purpose. In fig. |2| in the 
upper left plot the blue {green) region corresponds to the la allowed range for sin 2/3 {Rt) as 
calculated by means of (2.9). The red region corresponds to |e/^|^^ as given in (2.11). Finally 



the solid black line corresponds to a = 90° that is close to the one favoured by UT fits and 
the determination from B ^ pp |97j . 

The numerical input parameter that we use to obtain this plot are collected in tabs. [T] 
and [3] It is evident that there is a tension between various regions as there are three different 
values of (i?b,7), dependently which two constraints are simultaneously applied. 

Possible solutions to this tension can be achieved by assuming NP contributions just in 
one of the three observables ei^, sin 2/3, Rt. 
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Situation in the SM 



+24% NP effect in 




0.30 0.35 0.40 0.45 0.50 0.55 

-6.5° NP phase in Bd mixing 




0.30 0.35 0.40 0.45 0.50 0.55 




0.30 0.35 0.40 0.45 0.50 0.55 

-22% NP effect in AM<y/AM, 




0.30 0.35 0.40 0.45 0.50 0.55 



Rb 



Figure 2: The Rb — ^ plane assuming: i) sin 2/3, Rt and not affected by NP effects (upper left), 
ii) sin 2/3 and Rt NP free while €k affected by a positive NP effect at the level of +24% compared to 
the SM contribution (upper right), ill) ex and Rt NP free while sin 2/3 affected by a NP phase in 
mixing of —6.5° (lower left), iv) ex and sin 2/3 NP free while AM^/ AMg affected by a negative NP 
effect at the level of —22% compared to the SM contribution (lower right). The black star indicates 
the values for Ri, and 7 obtained in the NP UT fit of [83]. 



1) a positive NP effect in ex (at the level of +24% compared to its SM value) while taking 
sin 2/3 and AM^/ AMg SM-like [88], as shown by the upper right plot of fig. [2] The 
required effect in ex could be for instance achieved within models with CMFV by a 
positive shift in the function So{xt) [92] which, while not modifying {sm2(3)^Ks 
AMd/ AMs, would require the preferred values of sjEd^sFB^ ^ to be by ~ 10% lower 
than the present central values in order to fit AMd and AMs separately. Alternatively, 
new non-minimal sources of fiavour violation relevant only for the K system could solve 
the problem. 

2) ex and AMd/ AMg NP free while S^Xs affected by a NP phase in Bd mixing of —6.5° 



as indicated in (2.10) and shown by the lower left plot of fig. [2j The predicted value for 
sin 2/3 is now sin 2/3 = 0.823^q'q9q This value is significantly larger than the measured 
S-^Xs which allows to fit the experimental value of ex- 

3) ex and S^xs NP free while the determination of Rt through AMd/AMg affected by 



NP as indicated in (2.10) and shown by the lower right plot of fig. \2\ In that scenario 



one finds a very high value for \Vtd\ — 9.6 x 10 that corresponds to AM|^/AMf^ 



^ The value we give here is slightly smaller than the one found in |93l 1881 194) performing comprehensive 
fits of the UT. 
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P 


7] 


«[1 


7[1 


Rb 


iKfel X 103 


1) 


0.148t°:°i 


^29+°-°^^ 
u.oza_o.oi7 


93.2t^:° 


65.711° 


U.OUi_o 014 


44+0.17 


2) 




U.4Z4_Q 


86.6li^ 


65.711^ 


U.4UO_o.065 




3) 


U.UOO_o 062 


n Q79+0.023 
U.O(Z_o.oi9 


74.5t^r 


84.4l?o^, 


Q 074+0.021 
U-J '"±-0.018 


q p:7+0.23 
'-'•'J'-O.IQ 


UTfit 


0.177 ±0.044 


0.360 ±0.031 


92 ±7 


63 ±7 


0.404 ±0.025 


3.87 ±0.23 



Table 2: Predictions of several CKM parameters in the three scenarios as discussed in the text. For 
comparison, in the last hne results are also shown from a global NP fit of the UT [83] . 



(3.66^053) X 10~^, much higher than the actual measurement. In order to agree exactly 
with the experimental central value, one needs a NP contribution to AMj^/AMg at 
the level of —22%. Non-universal contributions suppressing AM^ {C'Ba < 1) and/or 



enhancing AM^ {Cb^ > 1) could be responsible for this shift as is evident from (2.10). 
The increased value of Rt that compensates the negative effect of NP in AMd/AMg 
allows to fit the experimental value of ex- 

Possibility 3) has not been discussed in [031 EHl El] • Interestingly, there are concrete and 
well motivated SUSY extensions of the SM compatible with scenario 3) as, for instance, abelian 
flavour models (see the following sections) or SUSY GUTs with right-handed neutrinos. In 
such cases, if the 6 — )■ s transition contains, in addition to a large mixing angle, also a natural 
C(l) CPV phase, then solution 3) also implies a non standard value for S^^. 

From fig. [2] it is clear that each of the solutions corresponds to particular values of Rb and 
7. In tab. |2]we show the values of the relevant CKM parameters corresponding to each case, 
where the values of the two variables characteristic for a given scenario are assumed not to 
be affected by NP. We observe 

• Solution 1) corresponds to 7 ~ 66°, Ri, ^ 0.36 and a ^ 93° in accordance with the 
usual UT analysis. 

• Solution 2) is characterized by a large value of R^ — 0.47, that is significantly larger 
than its exclusive determinations but still compatible with the inclusive determinations. 
The angles 7 ~ 66° and a ~ 87° agree with the usual UT analysis. 

• Solution 3) finally is characterized by a large value of 7 ~ 84° and a much below 90°. 
The latter fact could be problematic for this solution given the improving determinations 
of a. 

As seen in tab. |2l these three NP scenarios characterized by black points in fig. [2] will be 
clearly distinguished from each other once the values of 7 and Rb from tree level decays will 
be precisely known. Moreover, if the future measurements of {Rb, 7) will select a point in the 
Rb — 'j plane that differs from the black points in fig. |2| it is likely that NP will simultaneously 
enter ex, S^Ks a^id AMd/ AMg. 

We also note that the easiest way to solve the tensions in question is to require a particular 
sign in the NP contribution to a given observable: positive shift in |ex|, (t>Bd < and negative 
shift in AMd/AMg, with the latter implying increased values of Rt and 7. 

On the other hand, a negative NP contribution to ex would make the tensions in the 
Rb — 1 plane more pronounced, requiring stronger shifts in and Rt than in the examples 
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parameter 


value 


parameter 


value 


Bk 


0.724 ± 0.008 ± 0.028 [89] 


mAmt) 


(163.5 ±1.7) GeV [Ml EH] 


Fbs 


(245 ± 25) MeV HOOl 


mJmc) 


(1.279 ±0.013) GeV [101| 


Fb 


(200 ± 20) MeV fim 


Vcc 


1.44 ±0.35 [HSllTnS] 


Fk 


(156.1 ±0.8) MeV |103| 


Vtt 


0.57 ±0.01 [Bl] 

E 3 




1.22 ±0.12 llOOl 


Vet 


0.47 ±0.05 [Ml EH 1102] 


Be 


1.22 ±0.12 llOOl 

! i 


Vb 


0.55 ±0.01 [MII104] 


FbsVBbs 


(270 ± 30) MeV [MI] 


A 


0.2258 ± 0.0014 [H| 


FbVBb^ 


(225 ± 25) MeV [lOO] 


A 


0.808 ±0.014 [g 




1.21 ± 0.04 [100( 


Q 


0.177 ±0.044 [8] 


Vcb 


(41.2 ± 1.1) X 10"^ [S] 


V 


0.360 ±0.031 [8] 



Table 3: Input parameters used in the numerical analysis. 

given above. This specific pattern of the tension in the first plot in fig. |2] points towards 
certain NP scenarios and rules out specific regions of their parameter space. 



2.4 The Bs System 

Since the Bs system is central for our investigations, let us recall some known formulae. First 
the time-dependent mixing induced CP asymmetry 



r(5,(t) ^^0)±r(5,(t) ^^0) 



S^^ sin(AMst) 



(2.15) 



where the CP violation in the decay amplitude is set to zero. Next, the semileptonic asym- 
metry is given by 



_ T{Bs ^ l+X) - T{Bs ^ l-X) 
~ T{Bs ^ /+X) ± T{Bs l-X) 



Im 



r 



12 



Mf2 



(2.16) 



where represents the absorptive part of the Bs mixing amplitude. The theoretical pre- 
diction in the SM for the semileptonic asymmetry improved thanks to improvements in 
lattice studies of AB = 2 four-fermion operators |105] and to the NLO perturbative calcula- 
tions of the corresponding Wilson coefficients [106^ 1107] . 

Both asymmetries are very small in the SM where they turn out to be proportional to 
sin2|/3s| with Ps — —1°. The latter phase enters the CKM matrix element Vts 



Vts = -\Vts\e 



(2.17) 



As a consequence, both and S*^^ represent very promising grounds where to look for NP 
effects. 

In order to study NP effects in and S^^, let us recall possible parameterizations of 
the NP contributions entering the AF = 2 mixing amplitudes [83l 1108 



M- 



{Bs\H!s\Bs) = {Ml.f'^ + {M^_ 



12 



12; 



{{Ml^f^le^'f^^ + \{M- 



C 



■{M- 



\SM 



12; 



(2.18) 
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For the mass difference in the Bs meson system, one then has 



AM, = 2|M- 



121 



AM, 



SM 



In the case of the time-dependent CP asymmetry one finds [80] 

S^^ = - sin [Arg(Mi*2)] = sin(2|/3,| - 20^, 



(2.19) 



(2.20) 



where we took the CP parity of the final state equal to +1. 

Concerning A^^, we recall that, in the presence of NP, A^^ is correlated with S^^p according 
to the following expression jl08] . 



A 



SL 



Re 



12 



M 



12 



SM 



c 



1p<f> J 



(2.21) 



where we have neglected small contributions proportional to Im(rj2/^i2)^'^ used |106j 



s \SM\ 



(2.6 ± 1.0) X 10 



-3 



(2.22) 



Note that even a rather small value of S^rp ^0.1 would lead to an order of magnitude 
enhancement of A^-^ relative to its SM expectation, since within the SM is predicted to 
be of order 10"^ 

An alternative formula for the — model-independent correlation, pointed out 
recently in |109] . uses only measurable quantities 



A: 



AF. 



SL 



AM. 



1 



(2.23) 



and can be used once the data on AF. improves. In writing (2.23), we have assumed that 
ATs > 0, as obtained in the SM using lattice methods. 

Concerning the experimental situation, HFAG gives the following value for the semilep- 
tonic asymmetry [1], 

(^SL)exp = (-3.7 ± 9.4) X 10-3 . (2.24) 

In the case of S^^, there have been several analyses of the data from CDF [191 HIO] and 
DO [20] on the Bg — )• decay. Taking into account additional constraints coming e.g. from 
the flavour-specific Bg lifetime and the semileptonic asymmetry, tensions with the tiny SM 



-SM _ oi I ^ 0.036 at the level of (2 - 3)a are found [211 Ull ^- In the 
following we will use the results given by HFAG [Ij. For the Bg mixing phase, they quote 



prediction S^^ = sin2|/3s 



1/9.1-0 



Bs 



and give the following range at the 90% C.L 
This corresponds to 



(j)Bs e [0.10,0.68] U [0.89,1.47] . 
= OMtoil and S^^ E [0.20, 0.98] 



(2.25) 

(2.26) 
(2.27) 



respectively. 
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Figure 3: Left: The complex plane. The thick solid circles show the constraint on the Bg 
mixing amplitude from the measurement of AM^/AMg assuming no NP in AM^, while the thick 
dashed circles correspond to the constraint coming from AM^ alone. The thin lines correspond to 
NP contributions to with the indicated fixed phases 9s, see (2.18). The white regions correspond 
to the 90% C.L. range for 5^,^ in (2.27). The SM value corresponds to the black point. Right: Same 
as in the left panel, but now in the S^^-AMg/ AMf^ plane. 



In fig. |3l we show the experimental constraints in the complex Mfg plane and the S^^- 
AMg/ AMf^ plane, analogous to fig. [l] in the case. Interestingly enough, they show 
that the scenario with maximum CP violation, in which the phase of (M{*2)np is ~7r/2, is 
perfectly allowed and it would imply S^ip < 0.7 after imposing the constraint on -Ram = 
{AMd/AMs)/{AMl^/AMf^) given in tab. Q If we do not fix the phase of (Mf2)NP then 
all values for S^^ in the interval [—1, 1] are obviously possible, still being consistent with the 
constraint on Ram- 



3 AF = 0, 1, 2 Transitions in a General MSSM 

In this section, we discuss NP effects arising in general SUSY scenarios both in flavour vio- 
lating and flavour conserving processes. In the former case, we consider AF = 2 and AF = 1 
transitions both in the K and in the B systems, in the latter, AF = transitions such as the 
electron and the neutron EDMs de,n are discussed. 

Concerning the AF = 2 transitions, we present the SUSY contributions to the Bs^ and 
mixing amplitudes that enter the predictions of CP conserving quantities such as the 
mass differences AMg^^ and AMx and CP violating ones such as the time-dependent CP 
asymmetries S^Ks: S^^j, and ek- We also investigate the constraints coming from — _D° 
mixing. 

Regarding AF = 1 transitions, we give expressions for the CP asymmetries in — >■ 57, 
B — >■ (j){ri')Ks, the rare decays Bg^d — ^ f^^f^~ and B~^ — )■ t^u and discuss the three T-odd CP 
asymmetries in B ^ K*fj,~^fi~ . 
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In a general MSSM framework, there are various NP contributions to the FCNC processes 
that we consider. In particular, in the case of flavour changing processes in the down quark 
sector, one has contributions arising from one loop diagrams involving charged Higgs bosons 
and the top quark, charginos and up squarks, neutralinos and down squarks and also gluinos 
and down squarks. 

Gluino loops typically give the dominant contributions if flavour off-diagonal entries in 
the soft SUSY breaking terms are present. Such off-diagonal entries can be conveniently 
parameterized in terms of Mass Insertions. Throughout the paper, we will adopt the following 
convention for the Mis 

Ml) = diag(m^) + m^6d , = diag(m^) + m^6u , (3.1) 

where M.'q are the full 6x6 squark mass matrices in the super CKM (SCKM) basis in the 
convention of [113] (see also appendix [b] for the convention used) and iv? is an average squark 
mass. The Mis are then further decomposed according to the "chirality" of the squarks 

= ( ,Y ) ■ 

Our convention for the trilinear couplings Au and Ad follows also |113j and is such that 

^\5u%-i = -^t {At + ^l*/tp) , m'(5i^)33 = -m, {A, + ^i*tp) , (3.3) 
with mtAt = ^ (^«)33 and rribAb = ^ (^4^)33. 

In the remainder of this section, we will report expressions in the Mass Insertion Approx- 
imation (MIA) for the SUSY contributions entering in the processes we consider, since they 
are more suitable to understand the physical results 

In our numerical analysis instead, we work with mass eigenstates and we do not make use 
of the MIA as the latter cannot be applied when the flavour violating mixing angles are 0{1), 
as it is predicted by many flavour models. 

3.1 AF = 2 Processes 

We begin our presentation with AF = 2 transitions induced in the MSSM. To this end, let 



us briefly recall, as already stated for the Bg system in (2.18), that the meson-antimeson 
oscillations are described by the mixing amplitudes m[^^ = (MlT-L^^^lM), with M = K^, 
Bd,s- Within the MSSM, the effective Hamiltonian has the form 

5 3 

= ^^Q^ + E ^^^^ + ' (3.4) 

i=l i=l 



■^All the expressions we quote are relative to the MIA assuming a scenario with degenerate sfermion famihes. 
The phenomenological imphcations in flavour physics arising from a hierarchical sfermion scenario have been 



recently addressed in |114j . 
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with the operators Qi given, in the case of Bg mixing by 



Q4 = (s"Pin(^^^ij&^) , 

Q5 = (s"Pi6^)(s^Pij6") , 



(3.5) 



where Pr^l = |(1 ± 75) and a, /3 are colour indices. The operators 0i,2,3 are obtained from 
Qi,2,3 by the replacement L ^ R. In the SM only the operator Qi is generated, because of 
the {V — A) structure of the SM charged currents. On the other hand, within the MSSM, all 
operators typically arise. 

In a MFV MSSM, the NP effects to the mixing amplitudes are quite small both in the 
low tan /3 regime j56j and also at large tan /3 |115[ 1116] . 

In the former case, the largest contributions to the amplitude m[^^ arise from the chargino 
and charged Higgs effects by means of Ci and C3. All the other Wilson coefficients (WCs) 
involve couplings that are highly suppressed by light quark Yukawa couplings. Moreover, it 
turns out that the chargino box contributions to C3 are the only non-negligible contributions 
sensitive to flavour diagonal phases. In particular 



9l 



16n- 



*\2 



ml I 

:> n/rA -9/1 y^^^) 



+ 



m? 



(3.6) 



C. 



167r' 



*\2 
ts) ~, 



8M^ 



(3.7) 



where = \fi\'^/rh'^, tp = tan/3, the loop functions are such that /i(l) = —1/12, /3(1) = 1/20 
and e is the well known resummation factor arising from non-holomorphic (t^ enhanced) 
threshold corrections |117[ I118[ 1119] 1115] 1120] I121j . The dominant gluino contributions read 



2a, ^Mg 
o ~ o J y-^gj 



with Xg = M~/m'^ and the loop function / given in the appendix. One has /(I) 
that e ~ as/Sn for a degenerate SUSY spectrum. 



(3.8) 
1/2 such 



As anticipated, (3.6) clearly shows that C i is not sensitive to flavour conserving phases 



while C3 can be complex for complex fx At. In ( |3.7[ ) we have omitted contributions sensitive to 
the phases of the combinations /i^M| and fi'^AtM2, as they vanish in the limit of degenerate 
(left-handed) squark generations as a result of the super-GIM mechanism. 

Moreover, we note that ex can receive sizable effects only through Ci, as the contributions 
to ex from C3 are suppressed by ml /ml and thus safely negligible. We find that the Wilson 



coefficient Ci in (3.6) has the same sign as the SM contribution, thus we conclude that for 
flavour diagonal soft terms one has |ei^| > 



.SMl 
-K I 



Within a MFV framework at large tan /3, there are additional contributions to AP = 2 
transitions stemming from the neutral Higgs sector |122l I116[ I115[ 1123] 1124] I120j . However, 



^Analogous formulae hold for the B^l and systems, with the appropriate replacements of the quarks 
involved in the transition. 
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at least for /i > 0, these contributions turn out to be highly constrained by the experimental 
limits on BR(i?s — )■ /i+yU~), hence, they can be safely neglected. 

In the following, we therefore focus on the effects that arise in the presence of non-MFV 
structures in the soft SUSY breaking sector. In such a setup, the leading contributions to the 
K^, Bd and Bs mixing amplitudes come from gluino boxe^ The relevant Wilson coefficients 
in the MIA read IBljFl 



Cl - -^,m%2]g?{xg), (3.9) 

lit 

Cl ^ -^[(032]'^?r^(a:,), (3.10) 

Cl ^ [(5^^)32(032] ^/f^x,) , (3.11) 

fit 

Cl ^ -|| [(5i")32(5r)32] ^J'^x,) , (3.12) 

lit 



where Xg = M'^/m? and the analytic expressions for the loop functions g'^i \ g^l^ and g^^^ can 
be found in appendix A For the limiting case of degenerate masses we find g^^\l) = -1/216, 
g^4\l) = 23/180 and ^^(1) = -7/540. 

In the above expressions, we omitted the contributions arising from the LR and RL MI 



because, as we will see in sec. |3.2.1 they are tightly constrained by BR(6 — )■ 57). 

The same argument does not apply to the s ^ d and h ^ d transitions as there is no 
BR(i? — )■ Xs7) analog here However, in the concrete fiavour models we are dealing with, it 
turns out that the contributions from {Sd^)ij{6^^)ij are always dominant compared to those 
from and {Sd^)ij, hence, the above expressions are accurate enough for all the AF = 2 

transitions . 



An important observation comes directly from the Wilson coefficients (3.9)-(3.12): we 
expect large NP contributions in the mixing amplitudes from models that predict both nonzero 
iSd^)32 and (^1^^)32 Mis, since in this way the operators Q4 and Q5, that are strongly enhanced 
through QCD renormalization group effects [12711128] . have non- vanishing Wilson coefficients. 
We note in addition that especially the loop function g'^^ entering the Wilson coefficient C| 
is roughly a factor 30 larger than the one entering Cf. 



Expressions (3.9)-(3.12) are valid for the Bg mixing amplitude, while the corresponding 
expressions for B^ and mixing can be obtained by replacing the indices (32) with (31) 
or (21), respectively. Similarly, the corresponding expressions for — mixing are easily 
obtained from those relative to — making the replacement (5^)21 — ^ {8u)2i- 

In the case of mixing, it is important for our analysis to consider also the additional 
contribution coming from effective (2 — )■ 1) mass insertions generated by a double flavour flip 
(2 — )■ 3) X (3 — )■ 1). Results for the Wilson coefficients with one effective (2 — )■ 1) transition 

^In our numerical analysis we include the full set of contributions that can be found e.g. in ^96 . 
^The NLO Wilson coefficients in the MIA are given in [125] and in [1261 1121j in the mass eigenstate basis, 
^although, as a future perspective, LR and RL Mis relative to the 6 — >■ d transition will be probed by 
BR(B ^ p7). 

^Clearly the contributions of the Mis {S^^)ij and are taken into account in our model independent 

analysis of sec. |5] 
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are obtained in the third order of the MIA. We find 

Cl ^ -^,[^^',%M%M%i]9?\x,) . (3.13) 

Cl - -§ [(021(023(031] ^^^(x,), (3.14) 

Cl ^ m'M6r)2s{srhi + (5i")23(5i^)3i(5r)2i] (^.) , (s.is) 

fit Zi 

cl ^ m'M6f%s{6nsi + (5i^)23(5i^)3i(5r)2i] 9^'\x,) . (3.16) 

0(2) 
and for degenerate masses we find gl (1) = 

1/360, gf\l) = -1/6 and ^^^^(1) = 1/90. 

Finally we also have to consider the case where the (2 — )■ 1) flavour transition is entirely 
generated by (2 — t- 3) and (3 — ?■ 1) transitions. We find 

C! ^ -4[(023(03i]%fH^.), (3.17) 

^,2 



Cf - -^[(023(03l]^^r^(x,), (3.18) 

fit 

C'l ^ -||[(5i^)23(5,^^)3l(023(03l]^7f(a:,), (3.19) 

fit 

Cl ^ -^[(5i^)23(5i^)3i(023(03i]#(a:,), (3.20) 



with ^f^(l) = -1/3780, gf^{\) = 37/630 and g^^\l) = -1/378. The analytic expressions for 
these loop functions can again be found in appendix |Aj 

In the large tan /3 regime, the gluino box contributions are not the dominant ones anymore 
but they have to compete with double Higgs penguin contributions, see the Feynman diagrams 
of fig. |4| that are enhanced by tan^ /3 as in the MFV case |122l I116|, 1115] . Taking into account 
both gluino and chargino loops and generalizing the formulae of [122^ 1116^ 1115] to non-MFV 
contributions we find in the case of the Bg system]^ 

{C,)b - - 4^ (1 + etpY Mlm^ )32 [/^i(^.)J + 



2 2 +4 I |2 

Att 2M^ (1 + et,^)^ Mlm^ ^ 



(3.2i; 



with the mass ratios = |/ip/m^, X2 = \M2\^ /fh? and the loop functions, reported in 
appendix [A|, satisfying hi{l) = 4/9, /i3(l) = —1/4 and hi{l, 1) = 1/6. We stress that, due to 
the presence of the RR Mis, these contributions are neither suppressed by ms/mb nor by other 
suppression factors, as it happens in contrast in the MFV framework |122l 11161 11151 1124] . 

In the case of K'^ mixing, the most relevant effect from the neutral Higgses arises only at 
the fourth order in the MI expansion and we find the following expression 



'The Wilson coefficient for Bd mixing can be obtained by replacing (32) with (31) and Vts by Vtd- 
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sl bL sl 6i sl bi 



Figure 4: Feynman diagrams for the dominant Higgs mediated contributions to /S.B = 2 transitions. 
The leading contribution proportional to tan^ j3 comes from the self-energy corrections in diagrams 
where the Higgs propagators are attached to the external quark legs. 

with /i2 given in appendix [A| and /i2(l) = —2/9. 

In contrast to the case with gluino box contributions, now there are no analogous Higgs 
mediated contributions for — mixing; in fact, the non-holomorphic (tan /3 enhanced) 
threshold corrections lead to important FCNC couplings only among down-type fermions with 
the Higges. 



3.2 AF = 1 Processes 

3.2.1 Wilson Coefficients of the Dipole Operator 

The part of the AF = 1 effective Hamiltonian relevant for the 6 — )■ S7 transition that is most 
sensitive to NP effects reads 

•Heff D -^VttV* (C7O7 + CsOs + CrOy + CgOg) , (3.23) 

with the magnetic and chromomagnetic operators 

= ^m,{sa^^-PRh)F^, , Og = -|^m,(sa^^T^P^fe)G;^, . (3.24) 

Here, a^^ = |[7^,7'^] and T"^ are the SU{3)c generators. The operators Oi are obtained by 
the corresponding operators Oi by means of the replacement L ^ R. 

The dominant SUSY contributions to Cy^g arise from the one-loop charged Higgs, chargino 
and gluino amplitudes Cfg = Cj^^ + C^g + g. 

The charged Higgs contribution reads 

where yt = mf/Mjj±, e ~ 10~^ for a degenerate SUSY spectrum and the loop functions can 
be found in appendix |Xj 

The corresponding expressions for the chargino contributions to C7 and Cg read 



2 



r^X 

72 ~ m2 



Vt,V^*, m2 /7,8 1^2, X,) + ^2 ;7,g {X,) 



t 



/3 



[l + et 



/3J 



(3.26) 



The relevant loop functions are again defined in appendix \K\ 

Within a MFV SUSY framework, charged Higgs effects unambiguously increase the b ^ s'j 
branching ratio relative to the SM expectation, while the chargino ones (mostly from the 
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Higgsinos) can have either sign depending mainly on the sign and phase of fiAt. Since in 
SUGRA inspired models the sign of is set (in almost the entire SUSY parameter space) 
by the large (CP-conserving) RGE induced effects driven by the SU{3) interactions, the 
final MFV chargino effects will depend on the sign of the fi term that we assume to be real 
and positive which is preferred by the {g — 2)^ constraint (see also sec. 3.3.2). Then, the 



MFV contributions from the charged Higgs bosons and Higgsinos interfere destructively in 
BR(6 ^ 57). 

At the LO, both the charged Higgs and the chargino contributions to g are suppressed by 
ms/rrih. However, this mg/mf, suppression can be avoided in non-MFV scenarios at the NLO 
in the presence of RH currents by means of threshold corrections to the Yukawa interactions 
that are induced by tan /3-enhanced non-holomorphic effects |129] . For simplicity, we do 
not present these last contributions here, although they are systematically included in our 
numerical analysis. 

In the presence of non-CKM flavour structures, gluino mediated FCNC contributions also 
arise and the corresponding Wilson coefficients governing the 6 — )■ 57 transition read 



10 



4Gf 



'-'7,! 



9l 

TV? 



+ 



et,3) VtbV,l 



32 (2)/ 



(3.27) 



4G 



V2 



r^9 



frib VtbVt, 



32 (2)/ X 



+ 



{^d^)32 (3) 



(3.28) 



The functions ^f^' 



1,2) are given in appendix A As is evident from (3.27) and (3.28) 



;32 



and (0 



132 



b ^ s'-f puts very strong constraints on the helicity flipping mass insertions (0 
because the corresponding terms are chirally enhanced by a factor Mg/rrib- 

In fact, once the constraint from BR(6 — )■ 37) is imposed, these helicity flipping Mis cannot 
generate large effects in AB = 2 observables, in particular in S^^, anymore. 



3.2.2 Direct CP Asymmetry in b — >^ S7 

A very sensitive observable to NP CP violating effects is represented by the direct CP asym- 
metry in 6 — )■ 57, i.e. Acp{b — )• 57) [133] . If NP effects dominate over the tiny SM contribution 
A^{b — )■ 57) ~ —0.5%, the following expression for Acp{b — )■ 57) holds |134[ 1135] 



^ ' r(i?^x,7) + r(5^x,7) 



IC7 



(1.23 \m.[C2C;] - 9.52 ImfCgC;] + 0.10 ImfCaQ]) - 0-5 



m 7o 



(3.29) 



where we assumed a cut for the photon energy at ~ 1.8 GeV (see |134[ 1135] for details). 



In (3.29), the Wilson coefficients Ci are evaluated at the scale irib and they refer to the sum 
of SM plus NP contributions, i.e. C, = Cf^ + Cf^ (with Cf^ real). 

I n ord er to take into account effects from the Wilson coefficients Ci related to the operators 
Oj, (3.29 ) has to be modified according to CiCj — )■ CiCj+CiCj. Within a MFV scenario, Ci are 
completely negligible, being suppressed by a factor of rris/mb compared to the corresponding 



10 



Gluino contributions can also arise in the framework of the general MFV ansatz 0, see e.g. [1301113111132) . 



18 



Wilson coefficients Cj. This statement is no longer valid in the presence of NP in the right- 
handed currents, as is the case in some of the scenarios we are going to discuss. 

Still, as C*2 is negligibly small and the phases of Cj and are the same to a very good 
approximation in all the frameworks that we consider, Acp(& — )■ 37) receives NP contributions 
only through the imaginary parts of C^^. 

3.2.3 Time-Dependent CP Asymmetries in Bd — >■ cf){r]')Ks 

The time-dependent CP asymmetries in the decays of neutral B mesons into final CP eigen- 
states / can be written as 

Af{t) = Sf sin(AMt) - Cf cos(AMt) . (3.30) 

Within the SM, it is predicted with good accuracy that the 15*/! and Cf parameters are 
universal for all the transitions b — )■ q'q's {q' = c, s, d, u). In particular, the SM predicts that 
—TjfSf ~ sin 2/3 and ~ where rjf = ±1 is the CP eigenvalue of the final state /. NP 
effects can contribute tcE3 



i) the Bd mixing amplitude 

ii) the decay amplitudes b — )■ qqs (g = s, d,u) [8T], I136j . 

In case i), the NP contribution shifts all S/'s from sin 2/3 in a universal way while the C/s 
will still vanish. In case ii), the various S/s and also the C/'s are, in general, not the same 
as in the SM. 

The CP asymmetries Sf and Cf in B^ ^ f decays are calculated as follows. One defines 
a complex quantity A/, 

X^ = e-^'^''-^^^J(Af/Af) , (3.31) 

where is the NP phase of the Bd mixing amplitude, Mf2, and Af (Af) is the decay 
amplitude for Bd{Bd) — )■ /. Af and Af can be calculated from the effective Hamiltonian 
relevant for AB = 1 decays |137j . in the following way 

Af = {flHoslBd) , Af = {f\nos\B~d) , (3.32) 

where the Wilson coefficients of the effective Hamiltonian depend on the electroweak theory 
while the matrix elements {f\Oi\Bd{Bd)) can be evaluated, for instance, by means of QCD 
factorization fl38|. We then have 



2Im(A;) ^ _ 1 - |A; 



2 



The SM contribution to the decay amplitudes, related to 6 — )■ q'q's transitions, can always be 
written as a sum of two terms, Af^ = + A^, with A'j oc V;;Ks and oc IC^K.- Defining 
the ratio aj = e~^'^ {AJ / A'j) , we have 

Af^ = A} {1 + a}e'^) , (3.34) 

where the parameters have been evaluated in the QCD factorization approach at the leading 
order and to zeroth order in A/rrih in [138j . Within the SM, it turns out that S^fiKs — Sn'Ks — 



"^We assume that the asymmetry in the tree level transition h — ?> ccs is not significantly affected by NP. 



19 



S^Ks — sin 2(3 with precise predictions given in tabs. [T] and |6| The term provides only a 
neghgible contribution to Bd — )■ ipKs, thus A^^^ = — e"^*'^. Also for charmless modes, the 
effects induced by are small (at the percent level), being proportional to \(yubV*g) / (ycbV*g)\. 



The modification oi Af from the SM expression (3.34) due to NP contributions can always 
be written as follows 



12 



Af = A' 



1 + a]e'^ 



w, 



(3.35) 



where Cf^{Mw) and Cf^{Mw) are the NP contributions to the Wilson coefficients evaluated 
at the scale My/, the parameters and h'j^ calculated in |138j and C = ±1 depending on the 
parity of the final state; for instance C = 1 fo^' and C, = —1 for r]' Ks [139J. 
The generalization of this formalism to Bg decays is straightforward. 



3.2.4 CP Asymmetries m B ^ K*fi+fi 

The rare decay B — )■ K*{^ A'7r)yU+/i^ represents a very promising channel to look for NP in 
the B system, since its angular decay distribution is sensitive to the polarization of the K* 
and gives access to many observables probing NP effects. Furthermore, for the neutral B^ 
decay, the CP parity of the initial state can be unambiguously determined by measuring the 
charges of the kaon and pion in the final state. This "self-tagging" property allows a very 
clean access to CP-violating observables. 

In this paper, we focus on those CP-violating observables with relatively small dependence 
on hadronic quantities and large sensitivity to NP, following closely the analysis of [140] ; other 
recent analyses can be found in |14H I142j . In particular 

• all the observables are evaluated in the dilepton mass range 1 GeV^ < g2 < 6 GeV^; 

• power-suppressed corrections are accounted for by means of the full QCD form fac- 
tors in the naively factorized amplitude, and the "soft" form factors ^±^\\ in the QCD 
factorization corrections; 

• we focus on those observables which only depend on the ratios of form factors; 

• we assume NP effects in C/^g^i^s.p and C7^9,io,s,p but not in the other C,. 

In our analysis, we consider the T-odd CP asymmetries A-j, As and Ag which are not sup- 
pressed by small strong phases and can thus be sizable in the presence of new sources of CP 
violation (CPV) |141t 1140] . We use the conventions of |140j to which we refer for the full 
expressions of the asymmetries. 

In tab. |4| we just recall the main sensitivities of A-j, As and Ag to the relevant WCs of the 
AF = 1 effective Hamiltonian. Notice that (Ag) is only sensitive to Cy, Cg and Cio hence it 
represents, in principle, a golden channel to probe right-handed currents. However, as we will 
discuss in sec. [6| within the SUSY flavour models we are dealing with, the attained values for 
(Ag) seem to be far below the expected future experimental resolutions. This also implies that 
any potential signal of NP in (Ag) would disfavour the flavour models under study pointing 
towards different scenarios. 

^^We thank Dominik Scherer for pointing out the correct expression for Af |120j . 
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process 


CP asymmetry 


sensitivity to the WCs 


B K*fi+fi- 


Ar 


C7, C7, Cio, C*io 


B K*fi+fi- 


As 


Uj, 07, 09, 09, Olo 


B -> K*fi+fi- 


Ag 


C7, C9, Cio 


b — 7- 57 


Acpib -> S7) 


C7, Cs 


Bd ^ ^Ks 


S4,Ks 


Cs, Cs 


Bd ^ v'Ks 


Sr,'Ks 


Cs, Cs 



Table 4: CP asymmetries of various AF = 1 channels and the Wilson coefficients they are most 
sensitive to. 



We would also like to emphasize that (Ag) is peculiar from an experimental point of 
view since, in contrast to (A^) and (Ag), it can be obtained from a one-dimensional angular 
distribution and should thus be accessible at current B factories. In the conventions of |140] . 
the CP-averaged differential decay distribution in the angle reads 



1 

P + f 



^^^^ = ^ [(^2) + (^3) cos(20) + (Ag) sin(2</.)] , (3.36) 

where {S2) = —Fl has already been measured by BaBar |143j and Belle |144] and (S'3) is a 
sensitive probe of new CP conserving physics in right-handed currents |145[ I142[ 1140] . 

3.2.5 Bs n+^i' 

The SM prediction for the branching ratio of the decay Bg — )■ fJ'^fi'^ is (using the results 
in [80] with updated input parameters) 

BR(5, ^ /i+/i")sM = (3.60 ± 0.37) x 10"^ . (3.37) 

This value should be compared to the present 95% C.L. upper bound from CDF |146] 

BR{Bs ^ ^V")exp < 5.8 X 10-^ , (3.38) 

that still leaves a large room for NP contributions!^ In particular, the MSSM with large tan [3 
allows, in a natural way, large departures of BR{Bs — )■ fj.^fi^) from its SM expectation jl49[ 

The most relevant NP effects in the MSSM are encoded in the effective Hamiltonian 

HeS = —CsQs — CpQp — CsQs — CpQp , (3.39) 

with the scalar and pseudoscalar operators 

Qs = mt,{sPnb){ii) , Qp = mkisPpb) {i-f,i) , (3.40) 

as well as the corresponding operators Qs and Qp that are obtained by the exchange L i-^ R. 
For the corresponding Wilson coefficients in the MSSM, one has to a very good approximation 

Cp~-C5, Cpc::^Cs, (3.41) 



13 



An unofficial Tevatron combination with DO data |147] yields an upper bound of 4.5 x 10 ^ |148] . 
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Figure 5: Feynman diagrams for the dominant Higgs mediated contributions to Bg — t- 
The leading contribution to the decay amplitude proportional to tan^ /? comes from the self-energy 
corrections in diagrams where the Higgs propagators are attached to the external quark legs. The 
diagrams in the first row correspond to (3.42), the one in the second row to (3.43). 



with 



Cs 



ITLi 



Ml4M^{l + etpy{l + eetp) 
a2as f^e 



+ 



Ml 



l + et/3)2(l + e£t/3) m2 
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hUXg 



(3.42) 



- Ml4M^(l + et,)2(l + Qt,) rh- ^''^'^^'^ ' ^^'^^^ 

The loop functions hi, and hi already appeared in the discussion of the double Higgs 
penguin contributions to Bg mixing. This hints to correlations between NP effects in Bg 
mixing and Bg — )■ fj^~^fi~ identified within the MSSM with MFV at large tan /3 in |116l I115j . 

In fig. [s} we show the relevant SUSY Feynman diagrams for Bg in the presence 

of new flavour structures. 

The branching ratio for Bg can be expressed in the following way 



BR(i.,^.V)^^^;i-#(|B|^(l-4|-)+WM . (3.44) 



where the expressions A and B are given by the two linear combinations of the Wilson 
coefficients 

A = '^^Cff + ms, (Cp -Cp) , B = (Cs - Cs) , (3.45) 

where we assumed Cio free of NP, which is approximately true in all the scenarios that we 
consider. At leading order, the SM value for the Wilson coefficient Cio is given by 

C^iT = ^^^*.V;:Fo(xO, (3.46) 
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and the loop function Yq can be found in appendix |Xj The NLO QCD corrections to Yq have 
been calculated in |152j and found very small when mt{mt) has been used, which we do also 
in our analysis. 



3.2.6 T+u 

The SM expression for the branching ratio of the tree- level decay — t- r^z/ is given by 



BR{B+ T+u] 



SM 



Sir 



1 



mz 



m 



(3.47) 



Its numerical value suffers from sizable parametrical uncertainties induced by Fb+ and 146- 

On the theoretical side, the B ^ tu process is one of the cleanest probes of the large tan /3 
scenario due to its enhanced sensitivity to tree-level charged-Higgs exchange |153[ 1154^ 1155] . 
In particular, a scalar charged current induced by NP theories with extended Higgs sectors, 
leads to the following modification of the branching ratio 



R 



Btv 



BR(5+ ^ T+v) 
BR{B+ T+u) 



SM 



m 



B+ 



/2 



(l + et/3)(l + e^t/3). 



(3.48) 



where we have included the tan /3 enhanced non-holomorphic corrections for the quark and 
lepton Yukawas. In the limit of degenerate SUSY particles, it turns out that e ~ as/3n |115] 
and ei ~ — 3a2/167r |129] . 

Concerning the experimental situation, the HFAG collaboration [L\ quotes 



BR(5+ ^ r+z/; 



cxp 



;i.43±0.37) X 10" 



(3.49) 



which is based on results by BaBar |156j and Belle |157[I158] . Including additional preliminary 
results from BaBar |159j . one finds the following new World Average |112] 



BR(5+ ^ T+u] 



cxp 



;i.73±0.35) X 10" 



(3.50) 



that is considerably higher than (3.49). 

Using the value for |Kfe| quoted by the PDG |Kfe| = (3.95±0.35) x lO'^ and Fb+ ~ Fg 
given in tab. [3} we find the SM branching ratio given in tab. |6} corresponding to 



(^B..),.p = 1.57 ±0.53. 



(3.51) 



In view of the parametric uncertainties induced in (3.47) by Fb+ and Vub, in order to find the 
SM prediction for this branching ratio one can also use AMd to find 



BR(5+ ^ r+i/; 



SM 



4r]BSo{xt)BB,M^ 



1 - 



mz 



m 



B+ 



V, 



ub 



td 



tb+ AMd ■ (3.52) 



Here AM^ is supposed to be taken from experiment and \Vub/Vtdf is found using the formulae 

2Rt cos 13 



(ra, ^ and ^ 



Vub 


to 

1 


(1-AV2) 


Vtd 







(3.53) 
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with Rt and /3 determined by means of (2.9). In writing (3.52), we used Fb — Fb+ and 
~ mB+- We then find 

BR(5+ ^ r+z/)sM = (0.80 ± 0.12) x 10"^ (3.54) 



that is by roughly a factor of two below the data in (3.50). This result agrees well with a 
recent result presented by the UTfit collaboration |160j . 

It should be noted that the value of \Vub\ used effectively in this procedure turns out to be 
3.50 X 10~^, which is close to \Vub\ = (3.38±0.36) x 10~^ obtained from exclusive decays |161] . 
On the other hand, it is significantly lower than the one quoted by the PDG and obtained 
from tree level decays. 

For \Vub\ even higher than the PDG value, the experimental value for BR(i?+ — t- r+z/) can 
be reproduced, simultaneously making the appearance of a new phase in — likely as 



discussed in sec. 2.3[ In fact, the solution 2) to the UT tensions presented there would imply 



the central value BR(i?+ — )■ t~^i>) = 1.32 x 10~^ in the SM and be fully compatible with the 
experimental data, assuming no NP contributions to the latter decay. 

This discussion highlights the importance of accurate determinations of \Vub\ and of 
BR(i?"'' —7- r+z/) in the future in order to be able to decide whether NP is at work here 
or not. In particular, for low \Vub\ values (exclusive determination) there is a clear tension 
between the SM prediction for BR(i?+ — >■ t~^i>) and the data and an even larger tension exists 
in the presence of charged Higgs contributions. 

In our numerical analysis of NP contributions to various observables, we will use, to be 



conservative, (3.51) as the constraint coming from this decay. 



One could contemplate whether the charged Higgs correction in (3.48) could be so large 
that Rbtu becomes larger than unity improving the agreement of theory with the data. Such 
a possibility, that would necessarily imply a light charged Higgs and large tan (3 values, seems 
to be quite unlikely, even if not excluded yet in view of the constraints from other observ- 
ables [TUB] . 

This discussion is at first sight independent of the values of the weak decay constants that 



cancel in the ratio (3.52). Yet, the increase of the weak decay constants as suggested by recent 



non-quenched lattice calculations would certainly bring the SM value for BR(i?+ — )■ t^u) to 



agree better with the data. In (3.52), this effect is seen through the decreased value of \Vtd\ 
in order to agree with the experimental value of AM^ when the value of the relevant weak 
decay constant is increased. 

Another potentially interesting channel where to look for the presence of scalar charged 
currents is represented by the purely leptonic Kaon decays. In particular, the NP effect 



on RKfiu = T^^^^{K — )■ fii')/r^^{K — fiu) is obtained from (3.48) with the replacement 
— m|- |155] . Although the charged Higgs contributions are now suppressed by a factor 
m\-/m\ ~ 1/100, this is well compensated by the excellent experimental resolution |103j and 
the good theoretical control. The best strategy to fully exploit the NP sensitivity of 
systems is to consider the ratio R = Rkplv/ R■K^lv |155l 1103] instead of RK^lv■ In fact, while 
R and Rriiu have the same NP content (as Rttiiu is not sizably affected by NP), R depends 
on {Fk/Ft,Y instead of F|- with FkIF^, determined more precisely than Fx by unquenched 
calculations in lattice QCD. However, given that the resolution on Fx/ F^^ is at the % level, the 
same level of NP sensitivity of K — )■ Iv, we prefer to not include the constraints from K ^ iu 



^^For a detailed analysis on the viability of NP scenarios with a heavy-light extended Higgs sector, we refer 
the reader to Ref. |162j . 
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Process 


Present Bounds 


Expected Future Bounds 


Future Experiments 


BR(u ^6 7) 


1.2 X 10-11 


C(10-i3 - 10-14) 


MEG, PSI 


BR(u eee) 


1.1 X 10-12 


C(10-i3- 10-14) 


? 


BR(u ^ e in Nuclei (Ti)) 


1.1 X 10-12 


C(10-i8) 


J-PARC 


BR(r 6 7) 


1.1 X 10-^ 


c(io-8) 


SuperB 


BR(r -> eee) 


2.7 X 10-^ 


c(io-s) 


SuperB 


BR(r efifi) 


2. X 10-^ 


c(io-s) 


SuperB 


BR(r -> /i7) 


6.8 X 10-s 


c(io-s) 


SuperB 


BR(r fifi fi) 


2 X 10-^ 


o(io-s) 


LHCb 


BR(r /iee) 


2.4 X 10-^ 


c(io-8) 


SuperB 



Table 5: Present |81J and upcoming experimental limits on various leptonic processes at 90% C.L. 



in the present analysis. The above argument for K ^ £u does not apply to — )■ t^u. In fact, 
even if the hadronic uncertainties related to Fb and Vub are much larger that those for Fk/Ft^ 
and Vus, they cannot hide in any way the huge NP effects that can affect BR(i?+ — t- r+z/). 

3.2.7 K TTUU and b suu 

Within the MSSM with i?-parity conservation, sizable non-standard contributions to — >■ 
7rz/z7 decays can be generated only if the soft-breaking terms have a non-MFV structure. 
The leading amplitudes giving rise to large effects are induced by: i) chargino/up-squark 
loops |164l 11651 11661 1167] and ii) charged Higgs/top quark loops |168] . In the first case, large 
effects are generated if the trilinear couplings of the up-squarks have a non-MFV structure. 
In the second case, deviations from the SM are induced by non-MFV terms in the right-right 
down sector, provided tan/3 is large (30 to 50). 

In the case of 6 — )■ si'i' transitions like B — )■ Kuu, B — )■ K*^!) or i? — )■ X^z/z/, the second 
case above is prevented by the constraint on BR(i?s — > fi^fi^), while chargino/up-squark 
loops with non-MFV trilinear couplings in the up-squark sector can also generate sizable 
effects [T69lfT7n] . 

However, since the SUSY models we consider in sec. [6] feature neither sizeable off-diagonal 
entries in the trilinear couplings nor simultaneously large enough (5^^) 13 and {6^^)23 mass 
insertions, both K — )■ nuu and b — )• suu decays turn out to be SM-like. 

3.2.8 ii ijj 

Within SUSY models, LEV effects relevant to charged leptons originate from any misalign- 
ment between fermion and sfermion mass eigenstates. Once non-vanishing LEV entries in the 
slepton mass matrices are generated, irrespective of the underlying mechanism accounting for 
them, LEV rare decays like ii — ijj are naturally induced by one-loop diagrams with the 
exchange of gauginos and sleptons. The present and projected bounds on these processes are 
summarized in tab. IsIF^ 

^^The 2008 data from MEG are already close (BR(/i 67) < 3 x lO-" |171j ) to the present upper bound 
from MEGA so that the 2009 data should be able to provide a new improved bound. 
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The decay £i — )■ ij'y is described by the dipole operator and the corresponding amphtude 
reads 



(3.55) 



where p and q are momenta of the leptons ik and of the photon respectively and A^^ji are the 
two possible amplitudes entering the process. The lepton mass factor m^. is associated to the 
chirality flip present in this transition. The branching ratio of ii — )■ ij'j can be written as 



In the MI approximation it is found that |172j 



Al ^ 



Ol2 



47r m| 

2 



/iM2 



L(M|-/i2 



+ tan^ Ow yuMi 



fsnjXl) f2n{Xl,Xf,) 

m? ^ (^2 _ Ml) 



+ 



2 f2n{Xi) 



(3.56) 



4:71 



2 



(/x2 - Mf; 



+ 2 



Mi 
"^4 



f2n{Xi] 



, (3.57) 



where 6'vi/ is the weak mixing angle, mi is an average slepton mass, xi^2 = ^1,2/^% 
/x^/m? and fi(c,n){x, y) = fi(c,n){x) - fi{c,n){y)- The loop functions fi are given in appendix 



In the case of yU — )■ 67, one has to consider also the contributions arising from double 



'0 

Mis 



523531 as they can compete with the single MI contributions when ^21 ~ ^23531; in particular, 
this will turn out to be the case for the flavour models we will analyze. A particularly 
important effect is provided by the following amplitude 



m. 



^hn{x,){Sr)2,{Sf%, 



(3.58) 



because of the enhancement factor mT-/m^. The analog expression for Atj^ is obtained by 



A'i 



Af{L ^ R). In (3.58), the loop 



— 1/15 and flf\x) defined in appendix 



'unction is such that f^nix) = flf\x)/2 with /o'^''(l) = 
a| Other contributions, generated by the double Mis 
(5^''")23('5|'''')3i and {6^^)23(61^^)31 and not enhanced by rrir/rn^, can be still relevant; however, 
for simplicity, we do not report them here although they are systematically included in our 



<-{3), 



numerical analysis. In (3.56), (3.57), (3.58), as well as in the remainder of sec. we assume 



the /i term, the trilinear couplings and the gaugino masses to be real and the latter also 
positive. 

In the illustrative case of a degenerate SUSY spectrum with a common mass m^, we find 
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that 
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(3.59) 
(3.G0) 
(3.G1) 
(3.G2) 



Besides £j — )■ £^7, there are also other promising LFV channels, such as ii — )■ ijik^k and 
fi-e conversion in nuclei, that could be measured with the upcoming experimental sensitivi- 
ties. However, within SUSY models, these processes are typically dominated by the dipole 
transition ii — )■ ijY leading to the unambiguous prediction. 



BR{ii ijUjUi) ~ 37r 
CR(/i — )• e in N) ~ aorr 



log 



m 



4 



X BR(/i -> 67) . 



(3.G3) 



Thus, an experimental confirmation of the above relations would be crucial to prove the dipole 
nature of the LFV transitions. This would provide a powerful tool to discriminate between 
different NP scenarios as, for instance, SUSY and LHT models, as the latter do not predict 



a dipole dominance for 



[73]. 



Additional and sizable contributions to LFV decays may arise from the Higgs sector by 
means of the effective LFV Yukawa interactions induced by non-holomorphic terms |173] : 



hence, in general, the expectations of can be violated [T7il [T75l ITTGl [T771 fT78l fTTOllTSO] . 

However, these effects become relevant only if tan /3 = 0(40 — 50) and if the Higgs masses are 
roughly one order of magnitude lighter then the slepton masses |174l 11751 117GI 11771 11781 1179[ 
1180] . The last condition never occurs in our scenarios, hence Higgs mediated LFV effects are 
safely negligible in our analysis. 



3.3 AF = Processes 
3.3.1 Electric Dipole Moments 

As the SM predictions for the electric dipole moments are very far from the present exper- 
imental resolutions, the EDMs represent very clean probes of NP effects. Given that the 
EDMs are CP-violating but flavour conserving observables, they do not require in principle 
any source of flavour violation, hence, we refer to them as AF = processes. Yet, they can 
also be generated by two AF = 1 transitions, in which case one refers to "flavoured" EDMs. 

Indeed, within a MSSM framework with flavour violating soft terms, large and potentially 
visible effects in the "flavoured" EDMs are typically expected jlSlj 11291 1182] . In partic- 
ular, when NP sources of flavour violation generating 6 — )■ s transitions are assumed, the 
chromo-EDM (CEDM) and the EDM of the strange quark are unambiguously predicted. 
Unfortunately, an issue which is still unclear at present is the impact of the strange quark 
CEDM and EDM on the EDMs of physical systems like the neutron or heavy atoms like the 
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Thallium or the Mercury. The main source of uncertainty comes from the evaluation of the 
relevant hadronic matrix elements that should be ultimately evaluated by means of lattice 
QCD techniques. As a result, it is not possible at present to correlate or to constrain CPV 
processes in 5-physics by means of the EDMs of physical systems that are induced by the 
strange quark (C)EDM. Therefore, in our analysis, we only monitor the predictions for the 
strange quark (C)EDM generated by sizable CP violating effects in the flavour observables. 

If in the future there will be theoretical improvements enabling us to relate in a reliable 
way the strange quark (C)EDM to physical quantities, we could make use of additional 
observables, i.e. the EDMs of several systems, to test the NP theory that is at work, provided 
some NP signals in CPV Bg systems would appear. 

In the following, we report the relevant expressions for the "flavoured" EDMs including 
the dominant beyond-leading-order (BLO) effects |18H I129| 1182] . In fact, as shown in |18H 
1129^ I182j BLO effects dominate over the leading-order (LO) ones in a large region of the 
parameter space, hence, their inclusion in the evaluation of the hadronic EDMs is essential. 

Although our numerical results have been obtained including the full set of contributions, 
for simplicity, we report the dominant contributions to the hadronic EDMs. 

The dominant gluino/squark contribution to the down-quark (C)EDMs is 

1^. = -f^^'^l^f.^lm m%,i6r)u] . (3.64) 

e J - Air m"^ 1 + ec^ 

where the loop functions satisfy /^(l) = {4/135, 11/180}. 

Similarly, the corresponding prediction for the up-quark (C)EDMs is 

{t- = . (3.66) 

where k = 2,3, A^^^^^ = A^^t and f^{l) = {-8/135, 11/180}. 

The first effects to the (C)EDMs appear at the BLO [I82l [HH [I29] and they are well 
approximated by 

IT' "^^^J^, - "iG^^mUmI: 3(l + et,)2 )3.J fuAVt) , (3.66) 

where /^,±(1) = {-7/9, -2/3}. 

Charginos contribute to the EDMs already at the LO but the corresponding (C)EDMs 
are suppressed by the light quark masses rrid^. At the BLO, a new effect proportional to 
is generated by the charged-Higgsino/squark diagrams leading to |182[ I18H 1129] 

where = n"^ /w? and fH±{X) = {~5/18, —1/6}. The full expressions of the loop functions 
fg{xg), fH±{x) and //^±(a;) are listed in appendix [A} For equal SUSY masses and /x > 0, it 
turns out that eij = e/3 = a;s/97r. 

So far, we have presented the dominant contributions to the quark (C)EDMs assuming 
the presence of right-handed currents, hence of RR Mis. Since in the present work we are also 
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interested in models with purely left-handed currents, it is useful to show also the dominant 
(one loop induced) contributions to the (C)EDMs arising within this scenario. They read 



where gjj±{l) = {2/15, 1/10}, with its complete expressions given in appendix [A} Eq. (3.68) 



shows that, in this case, the (C)EDMs are suppressed by the external light quark masses, 
in contrast to the case where also RR Mis are non-vanishing. The analogous expression for 
up- type quarks, namely dif , is of order d^J^^ /du^ ~ [m^/(m^)f:^)] x [Atjj,/ jj,'^] xtj^ and thus safely 
negligible. 

Passing to the leptonic sector, the dominant contribution to the electron EDM arises from 
the one-loop exchange of binos/sleptons, and the corresponding EDM is given as 

+ {5i%k{5f'')u{5f%i + UuM] , (3.69) 



where fc, / = 2, 3, = -m,{A, + fitfs)/mj and /3n(l) = -1/12. 

One of the most peculiar features of the flavoured EDMs is that they might be proportional 
to the heaviest fermionic masses rrit, rrib and rrir instead of the lightest ones, as it happens 
in the case of flavour blind phases. This huge enhancement factor can bring the (C)EDMs 
close to the current and future experimental sensitivities, providing a splendid opportunity 
to probe the flavour structure of the MSSM. 

The main obstacles to fully exploit the NP sensitivity of the EDMs is that experimentally, 
one measures the EDMs of composite systems, as heavy atoms, molecules or the neutron 
EDM while the theoretical predictions are relative to the EDMs of constituent particles, 
i.e. quarks and leptons, thus a matching between quarks and leptons EDMs into physical 
EDMs is necessary and this induces sizable uncertainties related to QCD, nuclear and atomic 
interactions. 

The quark (C)EDMs and lepton EDMs can be obtained starting from the effective CP-odd 
Lagrangian 

Ces = - 5Z ■ ^h^i^i - ^ ^-^9si>i{G ■ a)-f5^/ji + dj (V'iV'i) (V'j^TsV'i) + • • • , 

i=u,d,s,e,fi i=u,d,s i,j 

(3.70) 

where the first and the second terms of (3.70 ) are the fermion EDMs and CEDMs, respectively, 



while the coefficients Cij are relative to the dimension-six CP-odd four-Fermi interaction 
operators. 

Among the various atomic and hadronic EDMs, a particularly important role is played 
by the thallium EDM (dTi) and the neutron EDM (dn). They can be estimated as |183[ I184[ 
CHHl [1861 [187] 

d^i = -585 4 - e 43 GeV , (3.71) 
dn = il± 0.5) [l.4 {dd - 0.25 rf„) + 1.1 e (c/^ + 0.5 <)] , (3.72) 

where Ci°^ is given by a combination of the coefficients Cij \186\ 1187] . However, when c?ti and 
dn are generated by flavour effects, as in our case, the contributions they receive from 
are always very suppressed hence, safely negligible |129] . 
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3.3.2 The Anomalous Magnetic Moment of the Muon 



The possibihty that the anomalous magnetic moment of the muon (we define = [g — 2)^/2), 
which has been measured very precisely in the last few years jl88j . provides a first hint 
of physics beyond the SM has been widely discussed in the literature. Despite substantial 
progress both on the experimental and on the theoretical sides, the situation is not completely 
clear yet (see |189] for an updated discussion). 

Most recent analyses based on e^e~ data converge towards a 3(j discrepancy in the 10~^ 
range |189j : 



Aa, 



(3± 1) X 10 



-9 



(3.73) 



The possibility that the present discrepancy may arise from errors in the determination of the 
hadronic leading-order contribution to Aa^ seems to be unlikely, as recently stressed in |191] . 

The main SUSY contribution to usually provided by the loop exchange of 

charginos and sneutrinos |192j . The supersymmetric contributions to are correctly re- 
produced by the following approximate expression 



MSSM 



47r mf^ 



-t- tan^ 6w 



L(M|-Ai2) V2 
yuMi 1 

-hn[Xl) + - 



-f2nix2,Xf,)-f2cix2,X^ 



m- 



2(M2-/i2 



(3.74) 



In the limit of degenerate SUSY masses one can easily find that 



„MSSM 
1 X 10-9 



1.5 



tan/3\ /300 GeVV 



10 



sgn /i . 



(3.75) 



The most relevant feature of (|3.75|) is that the sign of ajf^^^ is fixed by the sign of the /i term 



given M2 > 0) so that the solution /i > is strongly favoured. 



3.4 Correlations between AF = and AF = 1 Processes in the 
Leptonic Sector 

In the following, we discuss the implications of a potential evidence or improved upper bound 
of BR(/i — )• 67) at the expected sensitivities of MEG, namely at the level of BR(yU — )• ej) > 
10~^^ j204] . In particular, we will exploit the correlations among BR(yU — >■ 67), the leptonic 
electric dipole moments (EDMs) and the SUSY contributions to {g — 2)^ |205j . Finally, we 
discuss the prospects for the observation of LEV signals in r decays |205j . The corresponding 
numerical analysis in a concrete model will be performed in sec. [6] 

3.4.1 (g - 2)^ vs. BRiii ijj) 

An observation of /i — )■ 67 would provide an unambiguous evidence of NP but, unfortunately, 
not a direct test of the LEV source, as BR(yU — 67) depends also on other SUSY param- 
eters like the particle masses and tan p. While the latter parameters should be ultimately 
determined at the LHC/linear collider experiments, it would be desirable to access them by 

The most recent r-based estimate of the muon magnetic anomaly is found to be 1.9 standard deviations 
lower than the SM prediction |190j . coming closer to the e+e" value. 
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observable 


SM prediction 


exp. current 


exp. future 




~ 0.036 [81] 


o.8ii°i^ m 


~ 0.02 [193| 




sin2/3 + 0.02 ±0.01 [2] 


0.44 ±0.17 [I] 


(2 - 3)% [TM] 




sin2/3 + 0.01 ±0.01 [2J 


0.59 ±0.07 [I] 


(1 - 2)% [194] 


Acp{b 57) 


(-o.44;°ii) % m 


(-0.4 ±3.6)% [T] 


(0.4-0.5)% fTM] 




(3.4^°:^)10-3 [MI] 








(-2.6+°:^)10-3 |140| 








ro.i+°}uo-3 |i4o| 






\de\ (ecm) 


~ lO-^'' [196] 


< 1.6 X 10-^' [197] 


~ lO-'^i [196] 


\dn\ (ecm) 


~ 10-32 [196] 


< 2.9 X 10-26 [198] 


~ 10-28 [196] 


BR{Bs /i+/i-) 


(3.60 ±0.37)10-9 


< 5.8 X lO-'' [146] 


~ 10-^ [199] 


BR{Bd ^ /i+/i-) 


(1.08 ±0.11)10-1° 


< 1.8 X 10-s [Tie] 




BR{B X,7) 


(3.15 ±0.23)10-^ [200] 


(3.52 ±0.25)10-^ P 




BR{B -> Xsi+e-) 


(1.59 ±0.11)10-^ |201] 


(1.59 ±0.49)10-6 ^2021 1203] 




BR{B Tu) 


(1.10 ±0.29)10-^ 


(1.73 ±0.35)10-^ [112] 





Table 6: SM predictions and current /expected experimental sensitivities for the observables most 
relevant for our analysis. The branching ratio of i? — )• Xs^'^l~ refers to the low dilepton invariant 
mass region, g2^^_ G [1,6] GeV2. For the SM prediction of BR(i? tv), see also (3.54): BR(i? 
Tu) = (0.80 ±0.12) X 10-"^. 



exploiting the NP sensitivity of additional low energy observables. In particular, since both 
{g — 2)^ and BR(£j — )■ £^7) are governed by dipole transitions, the SUSY contributions to 
these observables are well correlated and their combined analysis provides a powerful tool to 
get access to the related LFV source. 

For a natural choice of the SUSY parameters, tp = 10 and a degenerate SUSY spectrum 
at m = 300 GeV, it turns out that Aa^^^"^ ^ 1.5 x 10-^ and the current observed anomaly 
can be easily explained. Assuming a degenerate SUSY spectrum, it is straightforward to find 



the correlation between Aa^^^"^ and the branching ratios of 



3-7 



BR{fi 67) 
BR(r ^ /i7) 



2 X 10-12 
8 X 10-^ 



■ Aaf ^ 2 

3 X 10-9 





3 X 10-9 



21 



10-4 



32 



10-2 



(3.76) 



where we have assumed that the Mis )ij provide the dominant contributions to BR{ii — )■ 



Eq. (3.76) tell us that, as long as the {g — 2)^ anomaly finds an explanation in SUSY 



theories, BR{£i — )■ £^7) are predicted once we specify the LFV sources. 

We emphasize that the extraordinary experimental sensitivities of the MEG experiment 
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Observable 


Exp. bound on 


(5f'^)i,- in RVV 


BR(/i 67) 


l(^[")21 


< 3 X 10-^ 


~ (0.3-1) X 10"^ 


BR(r ^ 67) 




< 6 X 10-2 


~ (2 - 6) X 10-3 


BR(r ^ /i7) 




< 4 X 10-2 


~ (0.3- 1) X 10-^ 



Table 7: Bounds on the effective LFV couplings from the current experimental bounds 

on the radiative LFV decays of r and /i leptons (see tab. p|) by setting Aa^^^"^ = 3 x 10-^. 
The expectations for the (5^^)jj's within the RVV2 model jSH 1206^ [22] are reported in the 
last column. The bound on (5f^)2i scales as [BR(/x — )■ e7)ea;p/1.2 x 10"^^]^/^, xhe scaling 
properties for the other flavour transitions are obtained analogously. 



looking for /i — )■ 67 offer a unique chance to obtain the first evidence for NP in low-energy 
flavour processes. Should this be the case, several leptonic observables related to /i — )■ 67 are 
also likely to show NP signals, i.e. the {g — 2)^, the electron EDM de but also other LFV 
processes like /x — )■ eee and fi — e conversion in nuclei. 

In order to get a more concrete idea of where we stand, in tab. [7| we report the bounds 
on the Mis arising from the current experimental bounds on BR(£j — )■ £^-7) imposing 

10 ^, corresponding to the central value of the {g — 2)^ anomaly. Moreover, 
in the last column of tab. [7| we also show the expectations for the Mis within a 

non-abelian SU{3) model that we will analyze in detail in following sections: the RVV2 
model [nil 1206^ [22] • Interestingly enough, the expected experimental reaches of MEG (for 
fi — )■ 67) and of a SuperB factory (for r — 7x7) would most likely probe the RVV model, 
provided we assume the explanation of the {g — 2)^ anomaly in terms of SUSY effects. 



3.4.2 Leptonic EDMs vs. BR{£i £jj) 

Within a SUSY framework, CP-violating sources are naturally induced by the soft SUSY 
breaking terms through i ) flavour blind F-terms |196[ 1207] and ii ) flavour dependent terms 
|181j . It seems quite likely that the two categories i) and ii) of CP violation are controlled by 
different physical mechanisms, thus, they may be distinguished and discussed independently. 

In the case i), the corresponding CP- violating phases generally lead to large electron and 
neutron EDMs as they arise already at the one-loop level. For example, when tp = 10 and 
mi = 300 GeV it turns out that 

4 ~ 6 X 10~25(sin Of, + 10~2 sin Oa) e cm , (3.77) 

while in the case ii) the leptonic EDMs, induced by flavour dependent phases (flavoured 
EDMs), read 

de ~ 10-22 X lm{{6f%^{Sf%,)ecm. (3.78) 

One of the most peculiar features disentangling the EDMs as induced by flavour blind or 
flavour dependent phases regards their ratios. In particular, 

^ = — - flavour blind phases, 

d^, 

de S,=2,3lm((^f^)ifc(^f^)fci) ^ ^ , r-^YQ^ 

— = — TTTT-s. — ; flavour dependent phases . (3.79) 

d^ Im((5f«)23(5[^)32) 
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In the case of flavour blind phases, the current bound dg < 1.7 x 10~^^ecm [8T] imphes that 
d^i ^ 3.5 X 10~^^ecm. On the contrary, in the presence of flavour dependent phases, the 
leptonic EDMs typically violate the naive scaling and values for d^ > 2 x 10^^^ e cm are still 
allowed. 

Moreover, when the EDMs are generated by flavour blind phases, they are typically un- 
related to flavour violating transitions On the contrary, the flavoured EDMs are closely 
related to LEV processes as they both arise from LEV effects and their correlated study would 
help to reconstruct the flavour structure responsible for LEV transitions. 

We recall that the predictions for the leptonic EDMs within a pure SUSY see-saw model 
are highly suppressed |208l 12091 1210] , at a level well below any future (realistic) experimental 
resolution. 

After imposing the current experimental bound on BR(/i — c-f), it turns out that de ~ 
10~^^ ecm, irrespective of the details of the heavy/light neutrino sectors. On the contrary, 
when the see-saw mechanism is embedded in a SUSY GUT scheme, as SU{5)YiN, de may 
naturally saturate its current experimental upper bound. As we will show in sec|6| also the 
RVV model naturally predicts large (observable) values for de- 

Hence, any experimental evidence for the leptonic EDMs at the upcoming experiments 
would point towards an underlying theory with either new sources of flavour blind phases or 
new CP and flavour violating structures beyond those predicted by a MEV scenario. In fact, 
as shown by (3.69), a crucial ingredient to generate non-vanishing EDMs is the presence of 
right-handed mixing angles and hence of right-handed currents. The latter are unavoidably 
generated in SUSY GUT scenarios as S'f/(5)RN or SU{10) through the CKM matrix and also 
in a broad class of abelian and non-abelian flavour models. 

A simultaneous evidence for the electronic EDM de and of /i — > 67 at the MEG experiment, 
could most likely suggest the presence of flavoured CP violating phases. 

Noteworthily enough, the synergy of low-energy experiments, as the leptonic EDMs and 
LEV processes (like /i — )■ ej), that are in principle unrelated, provides an important tool to 
unveil the anatomy of the soft SUSY sector. 



4 Soft SUSY Breaking and FCNC Phenomena 
4.1 Preliminaries 

The still elusive explanation of the pattern of SM fermion masses and mixing angles consti- 
tutes one of the main issues of the SM: the so-called "SM flavour problem". In addition, if 
nature is supersymmetric, the 'flavour problem" acquires a new aspect: whenever fermions 
and corresponding sfermions have mass matrices which are not diagonalized by the same ro- 
tation, new flavour mixings occur at the gaugino-fermion-sfermion vertices, generally leading 
to unacceptably large contributions to ECNC and/or CP- violating observables, of which K'^ 
mixing, and fi ^ cj are the most problematic. 

The most popular protection mechanisms to suppress such unwanted contributions are 

• Decoupling. The sfermion mass scale is taken to be very high. Still, such a scenario 
may be probed at the LHC through non-decoupling effects such as the super-oblique 
parameters |211j . 

-'^''A relevant exception is represented by the FBMSSM where it has been shown that there exist correlations 
among CP and flavour violating transitions in the _B-meson systems and the EDMs j68j . 
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• Degeneracy. The sfermion masses are degenerate to a large extent, leading to a strong 
GIM suppression. Such degeneracy could naturally arise from models with gauge- 
mediated supersymmetry breaking (GMSB) - or with some other flavour-blind mecha- 
nism of SUSY breaking mediation - if the mediation scale is low, ^® 

• Alignment. The quark and squark mass matrices are aligned, so that the flavour- 
changing gaugino-sfermion-fermion couplings are suppressed [23] . 

• MFV. Flavour violation is assumed to be entirely described by the CKM matrix even 
in theories beyond the SM [5]. 

We recall here that the MFV symmetry principle in itself does not forbid the presence 
of flavour blind CP violating sources [Ml EH [Ml EH EH [Ml EQl [212 IM, hence a MFV 
MSSM suffers, in general, from the same SUSY CP problem as the ordinary MSSM. Either 
an extra assumption or a mechanism accounting for a natural suppression of these CPV 
phases is desirable. The authors of [5J proposed the extreme situation where the SM Yukawa 
couplings are the only source of CPV. In contrast, in [71], such a strong assumption has been 
relaxed and the following generalized MFV ansatz has been proposed: the SUSY breaking 
mechanism is flavour blind and CP conserving and the breaking of CP only arises through the 
MFV compatible terms breaking the flavour blindness. That is, CP is preserved by the sector 
responsible for SUSY breaking, while it is broken in the flavour sector. While the generalized 
MFV ansatz still accounts for a natural solution of the SUSY CP problem, it also leads to 
peculiar and testable predictions in low energy CP violating processes [71] . 

As discussed in the previous section, in the general SUSY framework it is useful to pa- 
rameterize non-MFV interactions in terms of the Mis with {A, B) = (L, R) and 
{i,j = 1, 2, 3) on which the present data on FCNC processes put quite severe constraints [EH 
153] 1172] . Large departures from SM expectations are obviously still allowed in such a model- 
independent approach. 

While a model-independent analysis gives a global picture of still allowed deviations from 
the SM, its weakness lies in the fact that the suppression of FCNC processes is not achieved 
by some symmetry principle but basically by flne tuning the Mis so that existing experimental 
bounds are kept under control. Moreover, such an approach does not address the question 
of the hierarchies of the quark mixings in the CKM matrix nor the hierarchies of the quark 
masses. Analogous comments apply to the lepton sector. 

Much more ambitious in this respect are supersymmetric models containing flavour sym- 
metries that relate the structure of fermion and sfermion mass matrices. Such symmetries, 
while being at the origin of the pattern of fermion masses and mixings, can at the same time 
provide the sufficient suppression of FCNC and CP-violating phenomena by means of the 
degeneracy or alignment protection mechanisms discussed above. Moreover, as we will see 
in the context of our numerical analysis, SUSY flavour models imply certain characteristic 
patterns of flavour and CP violation that can be conflrmed or falsifled with the upcoming 
experimental sensitivities. 

Supersymmetric models with flavour symmetries have been considered extensively in the 
literature. They can be naturally divided into two broad classes depending on whether they 
are based on abelian or non-abelian flavour symmetries. They can be considered as gener- 
alizations of the Froggatt-Nielsen mechanism |216] : the flavour symmetry is spontaneously 

-'^^In GMSB models with a high messenger scale, gravity-mediated contributions cannot be neglected. The 
phenomenology of such gauge-gravity "hybrid" models is discussed in \212\ I213j . 
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broken by the vacuum expectation value of one or more "flavon" fields $j and the hierarchical 
patterns in the fermion mass matrices can then be generated by means of the suppression 
factors (($j)/M)", where M is the scale of the breaking of the flavour symmetry and the 
power n depends on the group charges of the fermions involved in the Yukawa couplings, 
generating the mass terms. 



4.2 Abelian Flavour Models 

There is a rich literature on models based on abelian flavour symmetries ^Hj [2H [23 ESI 123 
ESlESlEDlEIlESlESlElESlEniEZ]. While the simplest case, where a single ^7(1)^. group is 
employed, is disfavoured since it leads to unacceptably large contributions to FCNC processes 
as ek and AM/^ [231 [37], more successful models are realized through the abelian flavour group 
U{1)fi X U{1)f2 [23] • In this last case, the tight FCNC constraints are naturally accounted 
for thanks to a precise alignment of the down-quark and down-squark mass matrices [23] . 

On the other hand, the most prominent signature of this class of models are typically large 
effects in — Z)° mixing. In fact, abelian flavour symmetries do not impose any restriction 
on the mass splittings between squarks of different generations therefore they are expected to 
be non-degenerate with natural order one mass splittings. 

In particular, a mass splitting between the first two generations of left-handed squarks 
unavoidably implies a (1 <H- 2) flavour transition in the up-squark sector of order {8^^)21 ~ A. 

This can be easily understood by recalling that the SU{2)i gauge symmetry relates the 
left-left blocks of up and down squark matrices, i.e. (M^)^^ and (Mj)^^ respectively, in such 
a way that (M^)^^ = V*{M'^)^^V^ . In turn, the expansion of this relation at the first order 
in A implies that 

{Mlf^i = [V*{Mlf^V^]^^ ^ (Mj)^^ + A (m^ - m?) . (4.1) 

Thus, even for {Mj)2i = 0, which is approximately satisfied in alignment models, there are 
irreducible flavour violating terms in the up squark sector driven by the CKM as long as the 
left-handed squarks are splitted in mass. This is opposite to the case of non abelian flavour 
symmetries which we discuss now. 



4.3 Non-abelian Flavour Models 

In contrast to abelian models, where there is a lot of freedom in fixing the charges of the 
SM fermions under the flavour symmetry, non-abelian models [311 HOI [HI [121 [131 SH HSl 
Ii6| Wf\ [i8l [IHl [501 [SH [52] are quite predictive for fermion mass matrices once the pattern of 
symmetry breaking is specified. 

Moreover, non-abelian flavour models predict very small NP contributions for (1 ^ 2) 
flavour transitions, since, if the non-abelian symmetry was an exact symmetry, then at least 
the first two generations of squarks are degenerate, which would lead to a vanishing NP 



contribution to — and — mixings. However, as we will find in sec. 6^, the 
violation of the symmetry is strong enough to produce large effects in the observable ex- 

There are many candidates for the flavour symmetry group Gp, each having distinct sym- 
metry breaking patters. In general, Gp must be contained in the full global U{3)^ symmetry 
group of the SM in the limit of vanishing Yukawa couplings. In particular, a lot of attention 
is received by models with a U (2) symmetry jlTl HE] motivated by the large top mass, and 
also models with SU (3) symmetry [50] which are additionally able to naturally predict an 
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almost maximal atmospheric neutrino mixing angle 623 ~ 45° and to suggest a near maximal 
solar mixing angle 612 ~ 30°. 

4.4 Running Effects in Flavour Models 

Since the flavour models discussed in the previous subsections predict the pattern of off- 
diagonal elements of the squark mass matrices at the GUT scale while for the calculation of 
physical observables, their values at low energies are relevant, we now discuss the renormal- 
ization group (RG) evolution of these elements. In particular, crucial questions are whether 
these off-diagonal elements are strongly reduced or enhanced in the running, whether they 
mix with each other and whether they are generated if they vanish at the initial scale. In 
short, the question is whether the textures of the squark mass matrices predicted by the 
flavour models are RG stable. We will disregard off-diagonalities in the trilinear couplings, 
i.e. LR and RL Mis, as well as the slepton sector in the following discussion. 

Studies of these running effects in the context of the MSSM with Minimal Flavour Viola- 
tion have been performed in |217l 169] . 

A close inspection of the relevant RG equations (RGEs) |218] shows that for and 
m|), i.e. the RR Mis all mixing terms are suppressed by 1st or 2nd generation Yukawa 
couplings and can therefore be safely neglected. In fact, neglecting 1st and 2nd generation 
Yukawa couplings, the RGEs for the off-diagonal elements of read at the one-loop level 



167r2|(m?,),, 't.' 2 [y^] {nllUS^3 + ^.s) + ^Khlh , (4-2) 



IQTr'^^iniDij 't^ 2 [yl) (m^),,(5,3 + ^,3) + ^h^h^^ , (4.3) 



where t = log{fi/^o). As can be easily seen, {mfjj^)i2 are RG invariant in this approxima- 
tion; we have checked that this holds numerically to an excellent approximation even if light 
generation Yukawas and two- loop effects are taken into account. 

Concerning the remaining entries, we flnd that their values at low energies are well ap- 
proximated by 

K)^3 :^ 0.87 (m^)°3 , (ml)^^ ^ 0.82 {ml,)^ , (4.4) 
(m^)^3 ^ (1 - 0.10? - 0.05 1^) (m^)°3 , (4.5) 

(m^)23 ^ (1 - 0.10? - 0.05?) (111^)23 , (4.6) 

where we have defined t = tan/3/50, and where quantities with superscript on the right-hand 
side denote the values at the GUT scale predicted by the particular flavour model considered, 
while those on the left-hand side are meant to be evaluated at the weak scale. 

To summarize, the off-diagonal squark mass matrix elements in the RR sector are reduced 
by at most 15%, they do not mix among each other or with LL Mis, and they are not generated 
by the running once they vanish at the GUT scale. 

The situation in the LL sector is different; there, also mixing takes place, and the elements 
can be generated by RG effects even if they vanish at the GUT scale. Of course, both these 



"'^^The trilinear coupling matrices hu^d and the soft masses uIq jj ^ used in this section correspond to the 
conventions of Martin and Vaughn |218j and are related to the trilinear couplings and to the soft masses 
in the convention of |113j through the relations hu^d — "^ud ^^'^ ™Q = '^Q' d — i^u d)'^ ■ ^^^o 
appendix [B] for additional details. 
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effects are suppressed by combinations of CKM elements, since they would be absent if the 
CKM matrix were diagonal. The RG equation for the off-diagonal elements of niQ reads 



d 



16vr2^(m|)ij '= 2 {yd,iyd,j) (m^)^^ + [yl) (myij(5i3 + 5^-^) + y'l{m^Q)ik\kj + yl{m^Q)kjKk+ 



dt 



(4.7) 



where \ij = V^*Vtj and we have again neglected light generation Yukawas, except in the first 
term, which in the case of (ij) = (23) is only suppressed by ys/yb, but unsuppressed by CKM 
angles and can therefore be comparable in size to the remaining terms. 

Consequently we find that, numerically, the low-energy values of the (mQ)jj are well 
approximated by the following formulae. 



(my 



13 



23 



2 ^o 

12 



(0.91 - 0.05?) (my;'3 - Ami, - 0.09 [A12 (m^Jg + A23 (m^ 
(0.91 - 0.05?) (my°3 - Ami, " 0-09 Ui (mj) ^3 + A13 (111^°^ + 0.02? (1x1^)^3 



K)i2^K)i2-^^i2-0-09 



^13 (mQ)32 + A32 (my 



13 



where 



Ami 



Xij (0.33 ml + 0.89 Ml^ + 0.03 Al - 0.14 M1/2A0) 



(4.8) 

':>)2: 

(4.9) 
(4.10) 

(4.11) 



assuming CMSSM-like boundary conditions for the gaugino masses, trilinear couplings and 
the diagonal elements of sfermion mass matrices. 

As in the RR sector, the Yukawa-induced reduction of the elements, cf. the first terms 



in (4.7) and (4.8 )-(4.10 ), is only sizable in the (13,23) sectors. The terms in square brackets 



describe the mixing among the LL elements, while Amf^ describes the CKM-induced gener- 
ation of LL Mis, which takes place even in a completely flavour blind situation at the GUT 
scale, such as in the CMSSM. 

To summarize, the off-diagonal squark mass matrix elements in the LL sector mix among 
each other and they can be generated even if vanishing at the GUT scale; however, these 
effects are suppressed by combinations of CKM elements. Mixing of RR elements into LL 
elements only takes place in the 23-sector and is suppressed by a factor ysyt/yt- 

Finally, let us also remind that the attained values for the Mis 6ij are renormalization 
scale dependent as the diagonal elements are strongly affected by RGE effects. 



5 Strategy 

One of our main goals in the subsequent sections is to investigate the patterns of flavour 
violation in flavour models, putting particular emphasis on 6 — )■ s transitions. To this end, we 
select specific flavour models showing representative flavour structures in the soft masses. As 
we will see in sec. [6| in order to generate large NP effects in S^^, sizable - at least CKM-like - 
right-handed currents (driven by non-vanishing RR Mis) are unavoidable. The main reasons 
for this are twofold: 1) right-handed currents are less constrained than left-handed currents by 
low energy 6 — )■ s transitions, especially 6 — )• 57 and 2) since LL Mis (and hence left-handed 
currents) even if not present at the GUT scale, are RGE generated at the low scale via 
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the CKM and the top Yukawa couphng (see sec. 4.4), non- vanishing right-handed currents 
guarantee the presence of the large NP contributions provided by the left-right AF = 2 
operators (Q4 and in ( 3.5[ )) that are strongly enhanced by QCD RGE effects and by a 



large loop function. 

Therefore, we consider scenarios with 

i) large 0{1) RR mass insertions, 

ii) comparable LL and RR mass insertions that are CKM-like, 

iii) only CKM-like LL mass insertions. 

In sec. [6]we present their predictions for both CP violating and CP conserving effects occurring 
in 6 — )■ s transitions requiring at the same time that the models we consider satisfy the FCNC 
and CP violation bounds set by the experimental values of ex, AM^, AM^i, de, dn, etc. It 
will also be useful to compare the results of such an analysis with the results of the model- 
independent analysis discussed in sec. [2j 

To the best of our knowledge, the current analysis represents the most complete analysis 
present in the literature in this subject as for the number of processes considered and the 
inclusion of all relevant SUSY contributions. In particular, 

• We perform a full diagonalization of the sfermion mass matrices so that we do not make 
use of the MI approximation method. In fact, the latter method cannot be trusted when 
the flavour violating mixing angles are 0{1), as is predicted by many flavour models. 

• We systematically include the full set of one loop contributions to the FCNC processes, 
namely not only the gluino contributions, but also the charged Higgs, the chargino and 
the neutralino contributions. 

• We systematically include the beyond-leading-order threshold corrections arising from 
tan /9-enhanced non-holomorphic corrections in the presence of new sources of flavour 
and CP violation hence, accounting also for FCNC effects driven by the neutral Higgs 
sector. 

• All the above contributions have been systematically included for a very complete set 
of AF = 2, AF = 1 and AF = processes. 

We believe that the inclusion of all these contributions and observables is crucial in order 
to try to understand the pattern of deviations from the SM prediction that may be found in 
future measurements. 

In order to increase the transparency of our presentation, we will now outline the strategy 
for our numerics that will proceed in five steps. 



Step 1: Bounds on Mass Insertions 

In order to get an idea of the size of departures from the SM expectations that are still 
allowed in the supersymmetric framework, we will perform a model-independent analysis by 
calculating the allowed ranges for the mass insertions In this approach, as usual, only 

one MI of a given "chirality" AB and relative to a given family transition ij will be switched 
on at a time. Consequently, the bounds so obtained are valid barring accidental cancellations 
among different contributions. 
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Step 2: An Abelian Flavour Model 

In this step, we analyse an abelian flavour model by Agashe and Carone (AC) with a U{1) 
flavour symmetry [3B] which predicts large right-handed currents with order one CPV phases 
in the 6 — )■ s sector. 

As discussed previously, flavour models with a single U{1) are typically disfavoured by 
the ex and AM^ constraints [23l |37] . However, this is no longer true in the case of the AC 
model, where a high degree of quark-squark alignment is realized by means of a particular 
localization of fermions in extra dimensions, suppressing unwanted FCNC effects. 

In particular, the pattern of the relevant flavour off-diagonal Mis at the GUT scale, as 
function of the Cabibbo angle A, is given by 

0\ \ 

^ A2 , ~ ^ e'^^- , (5.12) 

A^ ^ / \0 e-'^^ * / 

(e)i2^A, (Oi2^AS (5.13) 

where we have suppressed the 0{1) coefficients which multiply the individual elements of the 
matrices and (pR is a free parameter. 

In our analysis, the correlations between various observables will play the crucial role. We 
will show that in the AC model the most interesting correlations involve the CP asymmetry 
Sjjj^ that can be much larger than in the SM so that this model can accommodate the recent 
data from CDF and DO. Indeed, in this context it is important to ask whether the desire 
to explain the latter data automatically implies other departures from SM expectations. We 
will therefore calculate BR(i?s — )■ fi~^fi~), S^Ks^ ^^^/i and as functions of S^^. 

The first correlation in our list is of particular interest as the upper bound on BR(i?s — 
/i+yU~) should soon be improved by CDF, DO and LHCb, possibly even finding first events 
for this decay. Next, the dependence of S^Ks 'S'^^ will tell us whether large values of 
the latter asymmetry are compatible with S^Ks significantly lower than S^Ks as signalled by 
BaBar pj. Similar comments apply to the {g — 2)^ anomaly Aa^. In addition, we investigate 
the ratios BR(i?s — > fi^fi^)/AMs and BR(i?s — ^ /i'*'/i~)/BR(-Bd — > fJ^^fi^) emphasizing the 
powerful tool they offer to unveil non-MFV structures of the model. We will also analyze 
the CP asymmetries in S — )■ K*^~^fi~ and the impact of the — mixing constraint on 
the correlations listed above. Finally, we will address the question how this model faces the 
current UT tension discussed in sec. 12.31 



Step 3: Non-Abelian SU{3) Models 

We next analyse a non-abelian model by Ross et al. (RVV) based on an SU{3) flavour 
symmetry |5T]. First, we recall that the observed structure in the Yukawa couplings does not 
flx uniquely the Kahler potential hence, the soft sector is not unambiguously determined. In 
the following, we analyse a particular case of the RVV model to which we refer to as RVV2 
model |62]. Similarly to the abelian AC model, also the RVV2 model predicts large right- 
handed currents. More explicitly, at the GUT scale, again suppressing the C(l) coefficients, 
the expressions for the flavour off-diagonal entries in the soft mass matrices in the SCKM 
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basis read l62l 



^3 ^iu}us y0.5^i{uj„a-x+l33) 




* — e —e 

^RR ^ I _^3g-ia;„. ^ -^0.5g-i(x-/33) | ^ (5x4) 

^^^0.5g-j(<^„s-2x+/33) ^^0.5gi{2x-/33) 

where the parameters £ ~ 0.05 and e ~ 0.15 are defined, after the fiavour symmetry breaking, 
in order to reproduce the observed values for fermion masses and mixing angles. Moreover, 
the phases Uus, X ^"^^ /^s are set, to a large extent, by the requirement of reproducing the 
CKM phase; in particular, it turns out that Uus ~ —A [62J and (x^/^a) ~ (20°, —20°) (or any 
other values obtained by adding 180° to each) [52] • Additionally, it is found that [S2] 

{^ll%2^\\ (Oi2^A^ (5.16) 

The trilinear couplings follow the same symmetries as the Yukawas. In the SCKM basis, 
after rephasing the fields, the trilinears lead to the following fiavour off-diagonal LR Mis |62] 

e^e-^"^- ^ £2 I ^mfc. (5.17) 

^3gi(t^„«+2/33-2x) ^2g2i(/33-x) 

Similarly to the fiavour model considered in step 2, S^^ can be large so that the analysis 
can be done along the lines of the one done for the AC model making the comparison of both 
models very transparent. In particular, as is already evident from the structure of the Mis, 
NP enters — Z)° mixing and in a profoundly different manner in these two models. 

As the RVV model is embedded in a 5*0(10) SUSY GUT model, correlations between 
fiavour violating processes in the lepton and quark sectors naturally occur making additional 
tests of this model possible. 

In particular, in the following, we list the fiavour off-diagonal soft breaking terms of the 
leptonic sector arising in the RVV2 model [02] 

-\e^ * eyl'^e-^^^-f^^^ \ , (5.18) 

_1^2^0.5gi(x-/33) ^^0.5gi(x-/33) ^ 

5^"^ ~ I -\e''e ^ ^^o.5g-.{2x-/33) | ^ (5.19) 

|^^^0.5g*{2x-/33) ^^0.5g^(2x-/33) ^ 

while the leptonic off-diagonal LR Mis have the following structure [62] 

7^ 



e 



:3 gi(2/33-2x) 3^2 g2i(/33-x) ^ / "^0 



6^/ - I ^ I l^m.. (5.20) 



In (5.14), (5.151, in order to avoid accidental cancellations among different phases, we have set to zero an 



extra CPV phase, ^2], that is not constrained by the requirement of reproducing a correct CKM matrix. 
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Finally, we consider a second example of non-abelian SU (3) flavour model analyzed by 
Antusch et al. |52] to which we refer to as AKM model. In the AKM model, in contrast to 
the RVV2 model, there is the freedom to suppress arbitrarily the flavour changing soft terms 
{^d^)iji hence we take the limit where = at the high scale. Therefore, our results have 

to be regarded as irreducible predictions of the AKM model, barring accidental cancellations 
among different contributions to physical observables. In the following, we report the flavour 
structure for the soft sector of the AKM model, in the SCKM basis, suppressing the C(l) 
coefficients 



where the above expressions are the same both for the down squark and the slepton sectors, 
modulo the unknown 0{1) coefficients. Concerning the up squark sector, the relevant flavour 
mixing angle generating — D^ mixing is (5^^)i2 — A^, hence, we conclude that the — D^ 
constraints are safely under control also in the AKM model. 

An interesting feature of the AKM model is the presence of a leading 0{1) CPV phase in 
the 23 RR sector but not in the 12 and 13 sectors. This will turn out to be crucial to generate 
CPV effects in the Bg mixing amplitude. Moreover, we notice that, while the RVV2 model 
predicts that {6^^)23 — e, the AKM model predicts that (5^^)23 — ^"^ with e ~ 0.15. Even if 
we expect that the CPV effects in the AKM model are smaller than in the RVV2 model, the 
indirect constraints (mainly coming from ex, b ^ s'j and AMg/ AMd), that are different in 
the two models, can play a crucial role to establish the allowed NP room for CPV effects in 
Bg systems. Hence, in order to distinguish between the RVV2 and the AKM model, by means 
of their footprints in low energy processes, a careful and comparative analysis is necessary. 

Step 4: Flavour Model with Purely Left-handed Currents 

We will next turn to the predictions for various low-energy processes induced by 6 — >■ s 
transitions within flavour models predicting pure, CKM-like, left-handed currents, i.e. Sj^^ 7^ 
and 6^^ = (or better Sf^ -C 5^^), with a CPV phase in the b —> s sectoi^^ 

This study will give us general predictions for a broad class of abelian [US] and non- 
abelian |12] flavour models not containing RH currents. Concerning the low energy predictions 
for the abelian case, we can draw clear cut conclusions immediately: the marriage of the tight 
constraints from — D mixing (typical of abelian flavour models) with the absence of RH 
currents (preventing large CP violating effects in AF = 2 observables in the down sector) 

allows non-standard and testable effects to occur only in CPV observables related to D'^ — D 
mixing. 

^^We will refer to this model as SLL model in the following. 




(5.21) 



The trilinear couplings lead to the following flavour off-diagonal LR Mis [521 




(5.22) 
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A non-abelian scenario has very different low energy predictions. In sec. [6} we will analyse 
a non-abelian model based on the following LL Mis 




(5.23) 

and all other Mis put to zero. Within that framework we will confirm our general statement 



of sec. 3.1 that S^^ does not receive large NP contributions in this context. As the NP effects 
in AF = 2 processes in such models are relatively small and a large asymmetry S^^ cannot 
be easily generated, this asymmetry is not an interesting variable in this scenario. 

In contrast, large non-standard effects can be generated in AF = 1 observables like S^Ksj 
Acp{b — 7- 57) and the CP asymmetries in i? — )■ K* . In particular, as we will see, S^Ks 
can attain values low enough to explain the related anomaly. We will therefore calculate 
this time Aa^, Acp(& — )■ 37), S^^iks ^-iid the CP asymmetries in 5 — )■ K*ii~^fi~ as functions 
of S(i,Ks- The predictions of models with purely left-handed currents will be systematically 
compared with those of somewhat similar models, i.e. the FBMSSM and the MFV MSSM, 
considered in step 5. 

Step 5: Comparison with the FBMSSM and the MFV MSSM 

Finally we comment on the predictions for CP violating effects arising within MFV MSSM 
models. This issue was recently addressed in [68] in the context of the so-called flavour blind 
MSSM (FBMSSM) [01 ESI EHl EZl ISH] , where the CKM matrix is the only source of flavour 
violation and where additional CP violating phases in the soft sector are allowed. 

In the FBMSSM, universal soft masses for different squark generations are assumed. Such 
a strong assumption gets somewhat relaxed in the framework of the general MFV ansatz [5] 
where the scalar soft masses receive additional corrections. The most general expressions for 
the low-energy soft-breaking terms compatible with the MFV principle and relevant for our 
analysis read [SI EH] 



2 ~ 2 
mn = m 



1 + feiYt + fosYjYd + (fesYlYrfYl Y„ + h.c. 



(5.24) 



where m sets the mass scale of the left-left soft mass, while hi are unknown, order one, 
numerical coefficients. The small departures from a complete flavour blindness of the soft 
terms generate additional FCNC contributions by means of gluino and squark loops. The 
latter effects were neglected in the context of the FBMSSM [HS] since, in principle, they can 
be very small if the parameters bi are small and/or if the gluino mass is significantly larger 
than the chargino/up-squark masses. In this respect, the contributions to FCNC processes 
discussed in [68] can be regarded as irreducible effects arising in MFV scenarios. 

The natural question that arises is whether the findings of [HE] , remain valid in the general 
MFV framework. Actually, within a general MFV framework, there are even new CPV effects 
to both AF = 2 and AF = 1 transitions, hence, potential departures from the FBMSSM 



predictions could in principle be expected. In particular in (5.24) the SM Yukawa couplings 
are not necessarily the only source of CPV. While bi and 62 must be real, the parameter 63 is 
generally allowed to be complex [69] . 



However, our analysis in sec. 6.5 will confirm the general finding of [HH] that within 
SUSY MFV scenarios, large NP contributions can only be expected in AF = 1 and AF = 
transitions. 
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We will also see that the FBMSSM and the MFV MSSM bear similarities to the models 
with left-handed CKM-like currents discussed in step 4. 



6 Numerical Analysis 

In this section, we present the numerical analysis of our study following the five steps described 
in the previous section. To this end, we list a number of constraints we impose throughout 
the analysis: 

i) the data on flavour physics observables (see tabs. [l]and[6]), 

ii) mass bounds from direct SUSY searches, 

iii) requirement of a neutral lightest SUSY particle, 

iv) requirement of correct electroweak symmetry breaking and vacuum stability, 

v) constraints from electroweak precision observables (EWPO) . 

We further assume a CMSSM spectrum where, in the case of the model-independent analysis 
(Step 1) we take at the EW scale only one non-vanishing MI at a time. Hence our bounds are 
valid barring accidental cancellations among amplitudes from different Mis. On the contrary, 
when we analyze abelian and non-abelian flavour models, we impose the flavour structures of 
the soft terms at the GUT scale where we assume the flavour models are defined. We then 
run the SUSY spectrum down to the EW scale by means of the MSSM RGEs at the 2-loop 



level |218j . consistently taking into account all the effects discussed in sec. 4.4 

The usual CMSSM contains (assuming vanishing flavour blind phases) five parameters: 
Ml/2, ^0, tan/3 and the sign of fi. However, we take a positive fi as is preferred by 
both the b ^ s'y and muon anomalous magnetic moment constraints. Consequently only four 
parameters are involved. 



6.1 Step 1: Bounds on Mass Insertions 

Starting with the model-independent analysis, we derive the allowed ranges for the Mis 
(with A,B = L,R and i ^ j = 1,2,3) under the constraints listed above. We scan the 
values of the CMSSM parameters in the following ranges: M1/2 < 200 GeV, mg < 300 GeV, 
IAqI < 3mo and 5 < tan/3 < 15, where the bound on Aq is set to avoid charge and/or colour 
breaking minima |219j . This choice for the input parameters would correspond to squark and 
gluino masses of order itIq, Mg < 600 GeV, both possibly observable at the LHC. 



6.1.1 1-2 Sector 

The measurements of AMx and ex are used to constrain the (5^'^) shown in fig. |6| AMk 
and eK constrain the real and imaginary parts of the product {8d^^2i {^d^)2V respectively. 
In the case of AM/f, given the uncertainty coming from the long-distance SM contribution, 
we use the range -AM^^p < (AMi^)|^ + (AM^)SUSY < aM^^p where (AM^)|^ refers to the 
short-distance SM contribution. The measurement of e'/e could put an additional bound on 
Im[((5^^)2i] that is effective in the case of the LR MI only. However, given the large hadronic 
uncertainties in the SM calculation of e'/e, to be conservative, we do not use this bound. 
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Figure 6: Bounds on the various Mis {6^^)21 (with A,B = L, R) as obtained by imposing the 
experimental constraints from tabs, [l] and [6j in particular ex and AMk- 



Notice that the bound on the RR MI is obtained in the presence of a radiatively induced LL 



MI {6j^^)2i oc Vf'jVts, see (4.10). In the Kaon sector, the product {5^^)21 {^d^)2i generates 
left-right operators that are enhanced by the QCD evolution, a large loop function and by the 
hadronic matrix element. Therefore, the bounds on RR Mis are more stringent than the ones 
on LL Mis. For the same reason also the constraints in the other cases ('5i^)2i = {^(1^)21^ 
{^d^)2i i^d^) 21 particularly strong. 

Similar arguments apply to the {3^^)^^ with {A,B) = {L,R) where we make use of the 
recent experimental measurement of D° — Z)° mixing. The corresponding bounds on (5^^) 
are based on |(M/^)susy| < 0.02 ps~^ |220] and shown in fig. [Tj 

Concerning the bounds on (5^^)21 (see the upper left plot of fig. [?]), one would naively 
expect that they are almost the same as the bounds found for (5^^)21! after the constraints 
from eK and AM^ are imposed. In fact, the SU{2)i gauge symmetry implies that {M'^)2i = 
[V*{M^)^^V^]^^, hence, {6^^)21 — {5^^)21 for degenerate left-handed squarks of the first 
two generations. However, as discussed in sec. |4| within scenarios predicting non degenerate 
left-handed squarks the above argument is no longer true, therefore {S!^^)2i 7^ {S^^)2i and 
the bounds on (5^^)21 and (5^^)21 apply independently. In fig. [7| we consider such a general 
scenario. Obviously, all the bounds on the other Mis {Su^)^^ with AB 7^ LL are unrelated 
to the corresponding bounds for the down sector as no symmetry principle is at work here. 

6.1.2 1-3 Sector 

The measurements of AMd and S.^Ks constrain the modulus and the phase of the Bd mixing 
amplitude, respectively. The bounds on the various {Sd^)r^^ are reported in fig. [sj The 
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Figure 7: Bounds on the various Mis (5„ )2i (with A,B = L, R) as obtained by imposing the 
experimental constraint from — mixing and the constraints from tabs, [l] and [g] 



constraints on {S^^)^^ and ((5^^)^^ are different because of the contributions of the left-right 
operator (arising by means of the RGE induced (5^^)31 oc VtbV*^) that is effective in the 

i^d^)3i ^^^^ only. The constraints in the cases {Sd^)r^^ = (^1^^)315 ('^1^)31 i^d^)3i^ 
are particularly strong due again to the large NP contributions provided by the left-right 
operators that are strongly enhanced by renormalization group effects and by a large loop 
function. 



6.1.3 2-3 Sector 

In this sector, we can exploit a large number of constraints. In particular, they arise from 
AMs and AB = 1 branching ratios such as 6 — )■ S7 and b — ?• si~^i~ . 

On the other hand, we do not impose the bounds from the time-dependent CP asymmetries 
in Bd — > (pKs and B^ v'^s and the direct CP asymmetry in 6 — )■ S7, given the still rather 
large uncertainties. In fig. [9] we show the allowed regions for the mass insertions and indicate 
in addition the resulting values for 5*^0 with different colors. The following comments are in 
order: 

• In the S strong constraint arises from BR(6 — )■ 57) as the related NP amplitude 

can interfere with the SM one. Regions in parameter space with a large real part of 
(^i"^)32 that would be allowed by 6 — )■ S7 are excluded by the constraint from BR{B — >■ 
Xsi^i~), while AMs does not provide any further constraint. We observe that in the 
considered framework the possible values for S^cf, are rather moderate and lie in the 
range —0.1 < S^^p < 0.1 even for quite large values of {S^'^)s2 < 0.1. 



45 





In the S^^ case, the situation is very different compared to the S^^ case as now BR(fe — )■ 
57) is not so much effective since the related NP amphtude (arising from right-handed 
currents) does not interfere with the SM one. Also BR(i? — )■ Xsi'^i") is not effective for 
the same reason. Now, AMg plays the main role in constraining {Sll^)s2- As in the other 
cases, RGE induced effects generate at the low energy an effective MI (5^^)32 oc VtbV^*. 
The product (5|'^)32('^d^'^)32 generates left-right operators that are enhanced both by the 
QCD evolution and by a large loop function. Therefore, the bounds on RR Mis from 
AMs are more stringent than the ones on LL Mis. We find that huge effects in S^tp in 
the entire range from —1 to 1 are possible in this scenario. 

In the case (5^'^) 32 = (^J^^)32, very strong constraints on the Mis arise from AMg due 
to the large contributions it receives from the left-right operator. Furthermore also 
BR(6 — )■ 57) provides additional constraints given the rather light SUSY spectrum we 
consider here. Still in the remaining parameter space large effects in the Bg mixing 
phase in the range —0.7 < S"^,^ < 0.7 are possible even for {6^^)32 = {Sd^)32 — 0.05. 

Finally, in the last two scenarios with (5^'^) 32 or (5^'^) 32 switched on, BR(6 — )■ 57) is the 
main constraint as the NP amplitude realizes the necessary chirality flip for the dipole 
6 — )■ 57 transition without involving the bottom mass insertion. In fact, the 6 — )■ 57 
constraint is so strong that S.^^ cannot depart significantly from the SM prediction in 
these cases. 



In fig. 10, we show the values attained by 5*^0 as a function of the modulus of different 



Mis, i.e. 1(5^^)32! (plot on the left), |(5f^)32| (plot in the middle) and \{5^%2i6^^)32\^^^ 
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Figure 9: Bounds on the various Mis {6^ )32 (with A,B = L, R) as obtained by imposing the 
experimental constraints from the b ^ s transitions hsted in tabs. [T] and |6j The different colors 
indicate the resulting values for 5^,^. 



(plot on the right), assuming a CMSSM spectrum with the following ranges for the input 
parameters: mg < 1 TeV, M1/2 < 1 TeV, \Aq\ < Sttiq, and 5 < tan/3 < 50. 

This choice of the parameter space is different from that employed before to get the bounds 
on the various 6s. In fact, in order to find the maximum allowed values for Sip^, it is crucial 
to consider the large tan /3 regime - in order to generate large Higgs mediated effects to 
S^(l, - and to enlarge the ranges for the SUSY mass scale - to relax the indirect constraints 
from observables (especially b — )■ 57) decouphng faster than S^^j, with respect to the SUSY 
mass scale. 

The plots of fig. 10 are complementary to those of fig. |9] and they show that the most 
natural scenario where it is possible to get large values for S^rp even for small, CKM-like 
mixing angles is the third scenario where simultaneously S^^ 7^ and 6^^ 7^ 0. 



6.2 Step 2: Abelian Model 

In this section, we analyze the abelian flavour model by Agashe and Carone [36], presented in 
sec. |5j stressing, in particular, the correlations among various observables. Here and in the 
analysis of the other flavour models, we assume a CMSSM spectrum with the following ranges 
for the input parameters: mo < 2 TeV, M1/2 < 1 TeV, \Aq\ < 3mo, and 5 < tan/3 < 55. 
Concerning the unknown 0{1) coefficients multiplying the off-diagonal entries in the soft 
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Figure 10: S^^ vs. |(<5^^)32|, \{S^%2\ and \{6l^%2{S^%2\^/^ respectively. All the points fulfill 
the indirect constraints from flavour physics. 



masses (5.12), they have been varied in the range ±[0.5,2] in our numerical analysis. 

As discussed in sec. |4| models with alignment naturally account for the tight FCNC 
constraints from ex and AMk by construction. Consequently in this model the ex anomaly 
discussed in sec. 2 can only be solved by modifying the value of Rt as also {sm2(3)^Ks is SM- 
like because of the negligibly small NP contributions to the h ^ d transition in this model. 
Yet we should remark that the latter solution implies a value for the angle a significantly 
below its best determined value (see fig. [2]) and it is to be seen whether such a solution 
will remain viable. On the other hand, in contrast to ex, large effects in — mixing 
are expected, provided the first two squark families are non-degenerate, as we expect from 
naturalness principles. 

In our numerical analysis, concerning the squark mass matrices at the GUT scale, we 
impose a large splitting between the 1st and 2nd squark generation masses such that rriu^ = 
2mc^ = 2mo, then the GUT scale MI is effectively {6^^)21 ~ A. However, at the low scale, 
where we evaluate the SUSY contributions to the physical observables, {6^^)21 is significantly 
smaller than A. In fact, there is a degeneracy mechanism triggered mainly by the flavour blind 
SU{3) interactions that restores a partial degeneracy between the 1st and 2nd generation 
squark masses |22H I222[ IBB] . In particular, while the off-diagonal entries in the squark 
mass matrices stay almost unaffected during the running from the GUT to the low scale 



(as discussed in sec. 4.4), the diagonal masses, in contrast, get strongly renormalized by the 
SU (3) interactions: their GUT scale values m| {Mqut 



Ami m| {Mgut 



ml become 

at the low scale m\^{Mw) — Am^ + QMy^ and mc^(Miy) ~ mg + QM^^. As a result, the 
GUT MI (5^^)21 ~ A is typically reduced by one order of magnitude at the low scale and 
the constraints from — mixing can be easily satisfied even for squark masses of a few 
hundred GeV. This is in contrast with the results of a low energy approach were (5^^)21 ~ A 
holds at the low scale implying a lower bound on the squark masses of around 2 TeV |220l 1223] . 



In fig. 11, we show the resulting predictions for several observables as functions of 



In all the plots, the black points satisfy the constraints of tab. |6] while the orange points 
correspond to negative NP effects in AM^/AMg at the level of 15% — 25% times the SM one, 
able to solve the UT tension. 



In the first row of fig. 11, we show, from left to right, the allowed values for Bg — )■ 



^(f, anu VS. b^tf,, respectively. In particular, fig. 11 shows that large 



vs. 5^0, S^Ks vs. 8^4, and vs. S.. 

values for S^^ in the range — 1 < S^^ ^ +1 are still allowed while being compatible with all 
the constraints such as — Z)° mixing, BR(6 — )■ S7), -Raa/ and ex- Interestingly, the first 
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Figure 11: Predictions for various observables vs. S^^ in the AC model. Orange points correspond 
to negative NP contributions to AM^/AM^ at the level of 15% — 25% times the SM one, able to 
solve the UT tension. The green points in the plot of BR(i?5 — )• vs. BR{B(i — )• and 

AMs/AMf^ vs. BR{Bs show the correlation of these observables in the MFV MSSM. 



plot shows that large effects in S^^p predict a lower bound on BR(i?s — )■ 

> 0.3 (the converse is obviously not true). 



) at the level of 
The correlation 



BR{Bs > 10-8 for 

between S^^i, and BR(_Bs — )■ /i^/i") signals that the double Higgs penguin contributions (3.21) 
are responsible for the large effects in S^cf, that we find in this model. 

Moreover, the correlation between S^cf, and S(f,Ks shown in the second plot indicates that 
both asymmetries can simultaneously depart significantly from the SM expectations although 
an explanation for the — )• (pKs anomaly would unambiguously imply negative values for 
S^fp, in contrast with the present data, in particular when the UT tensions are to be solved 
in this model. However, even in this case S^Ks can be at most suppressed down to 0.55. 

The shape of the correlation between S^^ and Sj^i^Tg, i.e. the fact that a positive (negative) 
S^(l, implies an enhancement (suppression) of S^f^Ks with respect to its SM prediction can also 
be understood analytically. In the considered model, the NP effects in S^Ks dominantly 
induced by the Wilson coefficient Cs given in (3.28). As we consider only the case of a real 
and positive /i parameter, the sign and phase o: 
to S^Ks is fixed by {5|^^)r^^■ Concerning S, 



Cs and therefore also of the NP contribution 
the dominant NP contribution to the mixing 



amplitude Mfg is induced by the double Higgs penguin contribution (3.21). In particular we 
find that in most parts of the parameter space the largest contribution comes from a double 
penguin with one gluino and one Higgsino loop (see the diagram in the middle of fig. lil). The 



corresponding analytical expression is stated in the second line of (3.21 ). As both the trilinear 
coupling At and the loop functions hi and have a fixed sign in almost the entire parameter 
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space considered by us, the sign and phase of the NP contribution to and hence to S^p^ 
is again determined by {Sd^)^2 correlation in the second plot of fig. emergesj^ 

As the CDF and DO data clearly favour S^^ to be positive, we conclude that froni the 
point of view of the AC model 

• Either < < 0.15, S^Ks ^ 0-60 and the UT tensions will go away, or 

• S^^ can be as large as 1, the UT tension is solved in this model but S(f,Ks will eventually 
turn out to be SM-like. 



The third plot confirms the model- independent correlation between S^^ and A^^ |108l 1109] 
discussed in sec. 2.4| and shows that the dileptonic asymmetry normalized to its SM value 



can lie in the range -100 < A|l/(A^l)sm < 100 for -1 < S^^ < 1. 

In the second row of the plot, from left to right, we show Aa^ vs. S^^, — )■ 
vs. Bg —7- /i''"yU~ and B^ — )■ vs. AMs/(AMs)sm5 respectively. Interestingly, there are 

many points accounting for large (non-standard) values for S^^ while providing a natural 
explanation of the {g — 2)^ anomaly. Moreover, the AC model can predict very striking 
deviations from the MFV SUSY expectations for 'BIi{Bd — )■ /i"'"/i")/BR(_Bs — > /i^/i~) and 
BR(i?s — > /i+/i^)/AM5/(AMs)sM- Hence, the above two ratios represent very clean and 
powerful observables to test the MFV hypothesis. In particular, as the AC model predicts 
new flavour structures only in the 6 — )■ s sector, it turns out that only Bs — )■ /i'^yU" can depart 
from the MFV predictions. More specifically, the ratio BR(i?d — )■ /i"'"/i~)/BR(i?s — )■ fi'^fi~) is 
dominantly below its MFV prediction and can be much smaller than the latter. 

Moreover, within a MFV framework, the current experimental constraints from BIi{Bs — >■ 
fj''^fj'~) already imply a very small NP room in AMs/{AMs)sy[, typically at a level smaller 
than 10%. In contrast, the AC model can predict quite sizable effects in AMs/{AMs)sm, as 



it is evident from fig. 11, which might saturate the Ram constraints. 



Further, even if not shown in fig. 11, the CP asymmetry in 6 — )■ 57 is SM-like as it is 
basically not sensitive to right-handed currents. 

In this model, NP contributions to K ^ ttuu decays are very strongly suppressed and the 
resulting branching ratios are SM-like. Consequently, if future experiments on these decays 
will show large departures from the SM expectations, the AC model will be ruled out. 

Concerning the CP asymmetries in — )■ K*fi'^fi~, we have found only small effects, at 
most at the percent level, after imposing all the indirect constraints, especially those from 

— D mixing. 

Finally, we point out that in the AC model, as well as in many other abelian flavour 
models, large effects for the neutron EDM are also expected. In fact, the peculiar flavour 
structure of the abelian flavour models - predicting large mixing angles in the 12 up-squark 



sector - leads to the following order of magnitude value for the up-quark EDM du (see (3.65)) 



du 
e 



|^ 2TeV y |lm [(eUQ^i] | ^ 



where we have assumed a degenerate SUSY spectrum, for simplicity. A similar expression 



holds for the up-quark CEDM too. Since dn ~ du (3.72), we conclude that 0{1) phases for 



22 



We note that if the dominant contr ibution to Mfj came from gluino boxes (3.11 ) or from double penguins 



with two gluino loops (first line of (3.21 )), with {Sj^'")r^2 induced radiatively through RG effects, the correlation 
between and S^Ks would have the opposite sign. 
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{^u^)i2 {4>ur) of abelian flavour models lead to a dn close to its current experimental upper 
bound. 

Analogously, CPV effects in — D mixing receive the dominant contributions by 
Im [((^^^) 12 (5^^)12] • Hence, large CPV effects in D'^ — D mixing would imply experimentally 
visible values for the neutron EDM by means of the up-quark EDM. 

Further, in the AC model, large CPV effects in Bs systems unambiguously imply a very 



large strange quark EDM. In particular, one can find that (see (3.64)) 



hence, the current experimental bounds on dn already tell us that either the strange quark 
contributions to dn have to be very small, with a proportionality coefficient smaller than 
10"^, or that 0{1) CPV phases for (5|'^)32 {4>dL) and/or (5f^^)32 {4>dR) are not allowed (unless 
0dL = ^d/?)- In the former case, large CPV effects in Bs systems are still allowed while they 
are excluded in the latter case (unless (pdL = 0di?)- In this respect, a reliable knowledge of 
the order of magnitude for the strange quark contributions to dn-, by means of lattice QCD 
techniques, would be of the utmost importance to probe or to falsify flavour models with large 
RH currents in the 2-3 sector embedded in a SUSY framework. 

In summary the striking predictions of the AC model in case the S^^ anomaly will be 
confirmed by more accurate data are: 

• The enhancement of BR(5s — /U^/W") above 10~^, 

• Aa^ ^ 10^^, thereby providing a natural explanation of the {g — 2)^ anomaly. 

• BR(5rf /i+/i-)/BR(B^ /^"^/i") and BR(5^ fj.+fi-)/AMs possibly very different 
from the MFV expectations with the first ratio dominantly smaller than its MFV value. 
Note however that BR(i?d — )> can reach values by a factor of 10 larger than in 
the SM. 



• Sff)Ks ~ S^Ks 

• Simultaneously UT tensions can be solved through the shift in AM^/AMs at the prize 
of a rather low a ~ 75°. 

• Large effects in — mixing but very small effects in — mixing. 



• Large values for the neutron EDM - very close to the current experimental bound - 
generated either by the up-quark (C)EDM or by the strange-quark (C)EDM. In the 
former case, visible CPV effects in — mixing are also expected while they are 
not necessarily implied in the latter case. A correlated study of several hadronic EDMs, 
with different sensitivity to the up-quark (C)EDM and the strange-quark (C)EDM, 
would provide a precious tool to unveil the peculiar source of CPV that is generating 
the neutron EDM. 
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6.3 Step 3: Non-abelian SU{3) Models 



We now turn to the analysis of another very interesting flavour model, the model by Ross and 
collaborators (RVV) [51J, or more specifically its particular version RVV2 recently discussed 
in [62] • It is based on the SU{3) flavour symmetry, as was presented already in sec. [sj As 
already stressed in previous sections, the symmetry properties of non-abelian flavour models 
naturally account for degenerate squarks of the first two generations (at least) solving thereby 
the flavour problem related to the experimental constraints from sk and AM^- As we will see, 
quite large - but still experimentally allowed effects in ex - can arise. Thus, in contrast to 
the AC model, the first solution to the UT tensions, through NP contribution to ex, becomes 
viable. Concerning the phenomenology of the up sector, the RVV model [51] and its RVV2 
version predict very small (most likely untestable) effects in the — mixing observables. 
These are probably the most peculiar features of non-abelian flavour models compared to 
the abelian ones with alignment, where the NP effects in ex are completely negligible by 
construction while large effects in — mixing are unavoidable. 

On the other hand, similarly to the AC model, also the RVV2 model exhibits large RH 
currents. Hence S^^f, can be large so that the analysis can be done along the lines of the one 
done for the AC model making the comparison of both models very transparent. 

A very interesting aspect of the RVV model (as well as its RVV2 version) is that it 
is embedded in a SUSY GUT SO(IO) model so that correlations among flavour violating 
processes in the lepton and quark sectors naturally occur making additional tests of this 
model possible. 

As a second non-abelian model, we also consider the SU{3) model proposed by Antusch 
et al. [52j (AKM model), already introduced in sec. |5] An interesting feature of the AKM 
model, that is in contrast to the RVV2 model, is the presence of a leading 0{1) CPV phase 
only in the 23 RR sector but not in the 12 and 13 sectors. This will turn out to be crucial to 
generate CPV effects in the Bg mixing amplitude. Moreover, in the RVV2 model it turns out 
that (6^^)23 — £, while in the AKM model we have {6^^)23 — £^ with e ~ 0.15. Hence, we 
expect larger CPV effects in the RVV2 than in the AKM models. We stress that, even if both 
models arise from an underlying S'f/(3) flavour symmetry, their low energy predictions can 
be quite different as their soft sectors, that cannot be set uniquely by the flavour symmetry, 
are different. 



6.3.1 RVV2 Model: Results in the Hadronic and Leptonic Sectors 



In fig. [T2| we show the predictions of the RVV2 model for several observables as functions of 
S^tl, in analogy to fig. 11 for the abelian flavour model. 



Our numerical results are obtained by implementing the flavour structures reported in 
sec. [5] assuming a CMSSM spectrum and scanning the unknown 0{1) coefficients in the range 
±[0.5,2]. 

Large values for S^^ up to —0.7 < S^^ ^0.7 are allowed while being compatible with all 
the constraints, in particular from BR(6 — )■ 57), -Ram and ex- 

One of the most prominent differences between the RVV2 and the AC models is the 
loss of correlation between S^^ and BR(i?s — )■ fi~^fi~), although both observables can differ 
spectacularly from their SM predictions. In particular, in the AC model, the Higgs mediated 
effects were the only contributions able to generate large (non-standard) values for S^tp, after 
imposing the indirect constraints, especially from b ^ s'j and — D mixing. As a result. 
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Figure 12: Predictions of the RVV2 model both in the hadronic and leptonic sectors. In the first two 
rows we show the predictions for various observables vs. S^^. The blue points correspond to positive 
NP effects in ex such that 1.2 < ex/(eA')sM < 1-3 and AM^^/AMs is SM-like. The green points in 
the plots oiBYl{Bs vs. BR(Bd Ai+/x") and AM,/AMf^ vs. BR(S^ show the 

correlation of these observables in the MFV MSSM. The last row refers to the predictions for leptonic 
observables. The green points explain the {g — 2)^ anomaly at the 95% C.L., i.e. Aa^j > 1 x 10^^. 
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Figure 13: Predictions of the AKM model both in the hadronic and leptonic sectors. In the first two 
rows we show the predictions for various observables vs. S^ff,. The blue points correspond to positive 
NP effects in ex such that 1.2 < ex / {^k)sm < 1-3 and AM^f/AMg is SM-like. The green points in 
the plots of BR(S^ ^ M^/^") vs. BR(5d ^ ^i+^x-) and AM^/AMf^ vs. BR(S^ Ai"^A*") show the 
correlation of these observables in the MFV MSSM. The last row refers to the predictions for leptonic 
observables. The green points explain the {g — 2)^ anomaly at the 95% C.L., i.e. Aa^ > 1 x 10~^. 
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S^p^ was correlated with Bs — t- fi^fi as the latter can also receive large SUSY effects only 
through the Higgs sector. 

In the RVV2 model, the absence of the D^—D constraints, as well as the more complicated 
flavour structure of the model (as for instance the presence of left-handed currents carrying 



new sources of CPV) destroy the above correlation as shown in the first plot of fig. 12 



The correlation between S^^ and S,f,Kg reveals that S^^f,, but probably not S'^^-g, can 
significantly depart from its SM expectation. Still, similarly to the AC model, negative NP 
values for S^Ks: &s suggested experimentally, imply negative values for S'^^, in contrast to the 
Tevatron results. 

Moreover, even if the RVV2 model predicts sizable left-handed currents with order one 
CPV phases, the CP asymmetry in 6 — )■ sj, as well as the asymmetries in i? — )■ K*i~^i~ turn 
out to be almost SM-like after imposing all the indirect constraints. 

In all the plots, all the points satisfy the constraints of tabs. [T] and [6] while the blue points 
correspond to a positive NP effect in ex such that 1.2 < eK/{€K)sM < 1-3 and AM^/AM^ 
SM-like, allowing to solve the UT tension discussed in sec. |2] However, it is interesting 
to observe that a simultaneous solution to the S^cf,, {g — 2)^ and ex anomahes is not very 
likely in this model, which constitutes an important test for this model. Similarly to the 
AC model the pattern of deviations from the SM predictions for Bg i^ — )■ /x+yU^ and AMg can 
differ spectacularly from MFV expectations. In particular, as in the AC model, the ratio 
BR.{Bii — )■ /i'*'yU~)/BR(i?5 — !■ yU'^/U") is dominantly below its MFV prediction and can be much 
smaller than the latter. 

On general grounds, we can conclude that the predictions of the AC and RVV2 models 
in the hadronic sector are quite similar. The NP effects in the RVV2 model are a bit smaller 
than in the AC model as the size of the right-handed couplings in the latter model is larger. 
However, one of the most important differences in predictions of these two models regards 
Ek and — mixing. We stress again that the RVV2 model naturally predicts sizable 
effects in ek (able to solve the current UT tension) while predicting negligibly small effects 
in — mixing; in the AC model exactly the opposite situation occurs. 

Concerning the predictions for the hadronic EDMs, we have explicitly checked that the 
down-quark (C)EDM reaches interesting values - up to ~ 10~^^e cm - while being compatible 
with all the constraints. Moreover, large CPV effects in Bg systems typically imply predictions 
for the neutron EDM within the expected future experimental resolutions ~ 10~^^e cm. The 
strange-quark (C)EDM is enhanced compared to the down-quark (C)EDM by a factor of 
^ e^^ ~ 50 hence, it would be very interesting to know precisely how it enters the neutron 
EDM. This would be also of great interest to establish the NP room left to CPV in Bg systems 
under the constraints from the hadronic EDMs. 

In the last row of fig. 12, we show the predictions of the RVV2 model for observables in 
the leptonic sector. From left to right of fig. 12, we report the correlations between S^^ vs. 
BR(/i — )■ 67), BR(/i — )■ 67) vs. the electron EDM and BR(r — )• fij) vs. BR(/i — )• 67), 
respectively. The green points are such that Aa^ > 1 x 10~^, thus they explain the {g — 2)^ 
anomaly at the 95% C.L. 

We observe that, both — )■ 67 and r ^ fi'j are very sensitive probes of LEV in the RVV2 
model and they both can turn out to be the best probes of LEV in this model. In particular, 
it is interesting to observe that the desire to solve the {g — 2)^ anomaly in the RVV2 model 
implies values of BR(r — )■ fi'-f) in the reach of LHCb and future SuperB facilities as well as 
BR(/i 67) > 10~^3, within the MEG resolution 



As we can see, both BR(£j — >■ ij'j) and de span over many orders of magnitude. Their 
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behavior can be understood looking at (3.56), (3.57) and (3.69), respectively. In fact, given 
that BR(£j — )■ ij'j) ~ [t^p/fh^) x \5ij\'^ and since tp and m can vary roughly by one order of 
magnitude in our setup while the |5jj|s are defined modulo unknown coefficients in the range 
±[0.5, 2], we expect BR(£j — )■ £^7) to vary by roughly seven orders of magnitude. 

Moreover, the loss of the correlation between BR(r — )• /i7) and BR(/i — )■ 67) can be traced 
back noting that BR(/i — >■ 67) receives several contributions (of comparable size) by both 821 
and ((52i)efr. ~ ^23(^31 while the only relevant Mis for BR(r — )■ /i7) are the ^32 ones. Further, 
the impact of the Mis and is much more relevant for BR(/i — )■ 67) - that receives 
an enhancement factor {m^/m^y' from the amplitude generated by and - than for 
BR(r — >• /X7) and this also contributes to destroy their correlation. 

Concerning de, we note that it is bounded from above by BR(/i — )■ 67) even if their 
correlation is loose given their very different sensitivity to the flavor structures of the model. 

In summary, the distinct patterns of flavour violation in the RVV2 model are 

• Large enhancements of S^^ and BR(i?5^ — )■ /i^/i^) but not in a correlated manner as 
in the abelian AC model, 

• Similarly to the AC model, BR(5d )/BR(5^ Ai^/W") and BR(5^ fi+fi-)/AM, 
might be very different from the MFV expectations with the first ratio dominantly 
smaller than its MFV value, especially for BK{Bs — ?► /i^yu^) ^ lO"*^. Also in this model, 
BR(i?d — 7- fi^fi^) can reach values by a factor of 10 larger than in the SM, 

• Removal of the UT tension through NP contributions to ex, requiring then typically 
S^p^ < 0.3 and S^Ks ~ S^Ksi 

• Small effects in — mixing, 

• Large CPV effects in Bg systems typically imply predictions for the neutron EDM within 
the expected future experimental resolutions dn ~ 10~^*e cm, 

• Large values for —0.25 < S^^ < 0.25 are still allowed even for BR(/i — )■ 67) < 10~^^. 
However, the desire of an explanation for the {g — 2)^ anomaly implies that BR(/i — >■ 
67) > 10-13 and l^^^l < 0.25, 

• BR(yU —7- 67) > 10^1'^, de > 10"^^ ecm and BR(r — )■ /i7) > 10"^ required by the solution 
of the {g — 2)^ anomaly. 

6.3.2 AKM Model: Results in the Hadronic and Leptonic Sectors 



In fig. 13, we show the predictions for the AKM model. As done for the other flavour 
models, we obtain the numerical results by implementing the flavour structures reported in 
sec. [5] assuming a CMSSM spectrum and scanning the unknown 0{1) coefficients in the range 
±[0.5,2]. 

The main differences between the RVV2 and the AKM models can be traced back re- 
membering the peculiar flavour structures in the soft sector of the two models. In particular, 
the AKM model predicts a CKM-like RH current while the corresponding mixing angle in 
the RVV2 model for the 6 — ?■ s transition is larger. This implies that in the AKM model, 
the effects in CPV observables are typically smaller, but still very interesting, compared to 
the RVV2 model. In particular, it is found that, in the AKM model S^prj, lies most likely in 
the range —0.3 < S^^ < 0.3 while being compatible with all the constraints. Interestingly 
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enough, similarly to the abelian case discussed before, large non-standard values for S^^ would 
unambiguously point towards non-standard values for BR(i?s — > as shown in the first 



plot of fig. 13 However, fig. 13 also shows that within this model the S^j^g anomaly cannot 
be accounted for. Large values for S^^ can also be compatible with an explanation of the 
{g — 2)^ anomaly. Moreover, we also observe that departures from the MFV SUSY expecta- 
tions for BR{Bd /i+/i")/BR(S, and BR(S, /i+/i-)/AM,/(AM,)sM(in both 
directions with respect to the SM predictions) are also expected in the AKM model, even if 
the ratio BR(i?rf — )■ /i+yU^)/BR(_Bs — )• /i^/i^) stays much closer to the MFV value of roughly 

1/33 pMig]. 

Passing to the hadronic EDMs, we observe that, even if the AKM does not contain leading 
CPV phases in the 13 sector, the down-quark (C)EDM is always generated by means of the 



CKM phase, in the presence of RH Mis (see sec. 3.3.1). We have explicitly checked that the 
down-quark (C)EDM might reach values up to ~ 10~^®e cm after imposing all the constraints. 
As a result, the constraints from the hadronic EDMs are well under control. However, large 
CPV effects in systems would likely imply predictions for the hadronic EDMs within 
the expected future experimental resolutions ~ 10~^^e cm. The strange-quark (C)EDM is 
enhanced compared to the down-quark (C)EDM by a factor of ~ ~ 50, as in the RVV2 
model. Once again, we stress that it would be of crucial importance to know how the strange- 
quark (C)EDM enters the hadronic EDMs in order to establish which are the CPV signals 
we can still expect in Bs systems. 



In the last row of fig. [T3| we show the predictions of the AKM model for observables in 
the leptonic sector as done for the RVV2 model in fig. [T2j In particular, both — )■ 67 and 
r — 7- /i7 are sensitive probes of LEV in the AKM model. However, in contrast to the RVV2 
model, /i — )■ 67 will represent the best probe of LEV in this model, especially after the MEG 
sensitivity will be fully exploited. In fact, it turns out that when BR(/i — )■ 67) ~ 10~^^ then 
BR(r — )■ /Z7) < 10"^*^, far from the SuperB reach. 

Moreover, the predictions for the electron EDM in the AKM model are well below those 
of the RVV2 model. In fact, in the RVV2 model, the two largest contributions to are 
proportional to either lv[i{5^^)i2,{5f^)^i or Im(5^^)i3(5|*^-^)3i with {5f^)^i carrying a lead- 
ing 0{1) CPV phase. In contrast, in the AKM model, the first non vanishing contribution 
to de is generated at the third order in the MI expansion by means of the combination 
Im(5|'^)i3(5f^)32(5f^)2i as (5f^)32 is the only MI containing an 0(1) CPV phase. Hence, 
there is a higher order suppression in terms of small mixing angles in the AKM model com- 
pared to the RVV2 model. 

An order of magnitude for the upper bound on rfg, compatible with the constraints from 
BR(/i — )■ 67), can be obtained considering a degenerate SUSY spectrum and assuming that 
the dominant contributions to BR(/i — )• 67) come from {5f^)2i and {5^^)2i- 



Taking the specific expressions for the Mis arising in the AKM model (5.21) and (5.22) 
one can easily find that 



<r in-29. /BR^^ 67) 
10 



~ \l Tc^-u sin ^ 6 cm , (6.3) 



as is fully confirmed numerically by fig. 13 



Still, visible values for rfg can be reached in the AKM model; in particular, it turns out 
that 4 < 10~29(iQ-30) _^ ^ 10-^(10-13). 

Finally, we note that, for BR(/i — )• 67) < lO'^^, only values for S^^ at the level of 
|>S'^(/)| < 0.15 are still possible if one wants to be fully compatible with an explanation for the 
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{g — 2)^ anomaly. 

Therefore, the disentanghng of the RVV2 and the AKM models might be problematic 
from a low-energy point of view, even if not impossible. Clearly the knowledge of some SUSY 
parameter, as the value for tan /3 and the charged Higgs mass, for instance, would be of 
outmost importance to make access to the flavour structure of the flavour model at work. 

In summary, several of the predictions in the AKM model are similar to the ones found 
in the RVV2 model, but the following most significant differences should be noted: 

• \S^(p\ can reach values up to |5'^</,| ^0.3, values that are fully consistent with an expla- 
nation of the UT tension through NP contributions to ex 

• S^ci, > 0.2 uniquely implies BR(_Bs — )■ > 10~^. 

• BR(i?rf —7- /i+yU^)/BR(i?s — 7- fi^fJ^") departs less from MFV expectation than in the AC 
and RVV2 models and can be both larger and smaller relative to its MFV value. 

• A non standard S^^ typically imphes dn ~ 10~^^e cm that is within the expected future 
experimental resolutions. 

• BR(r — fi'-f) is likely out of the reach of SuperB machines, in particular, BR(r — )■ 
Ai7) ^ (10-^ 10-^ 10-1°) when BR(/i 67) < ilO-^\ IQ-^^, 10~^^), respectively. 

• 4 < 10-29(10"=^°) ecm for BR(/i ^ 67) < 10-"(10-i3) ^j^g^e BR(/i ^ 67) > lO^^^ jg 
required by the solution of the {g — 2)^ anomaly. 



6.4 Step 4: Flavour Model with Purely Left-handed Currents 



Another typical flavour structure for the soft sector emerging in many abelian and non-abelian 
flavour models are purely left-handed currents with CKM-like mixing angles 



14 



Based on the flavour structures given in (5.23) and scanning the unknown 0{1) coefficients 
in the range ±[0.5, 2], we present in fig. 



the numerical results of our study of the non-abelian 

model of |12]. 

As already discussed in sec. [51 in the abelian case, non-standard and testable NP effects 
can be generated only in observables related to — D mixing hence, in the following, we 
focus on the non-abelian case. 

The first upper plot refers to the correlation between v4cp(& — )■ 57) vs. S^Ks- We find that 
—0.2 < S^Ks ^ 1 and —6% < Ac^ip — )■ 57) < 6%. Interestingly, it turns out unambiguously 
that positive (negative) NP contributions for S^Ks a-^e associated with negative (positive) NP 
effects in Acp(6 — )■ 57). 

This correlation can also be understood analytically. While the NP contributions to S^Ks 
arise dominantly from the Wilson coefficient Cg, the CP asymmetry Acp(& — )■ 57) crucially 



depends on the relative size of the imaginary p arts o f C7 and Cg, see (3.29). In the considered 
model, gluino loops typically give Cj < C| (3.27), while Wino loops lead to the opposite 



situation, i.e. C7 > Cg (3.26). In fact the correlation that we find corresponds to the 
latter case, implying that Wino contributions dominate over gluino ones in large parts of the 
parameter space, since Winos are typically a factor 3 lighter than gluinos. 

In the second upper plot we show that Sr^iKs turns out to be highly correlated with S^Ks-i 
as they both receive the dominant NP contribution (that is larger in the S^Ks case) from the 
chromomagnetic Wilson coefficient Cg. Moreover, a solution of the S^Ks anomaly can also 
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Figure 14: Predictions for various low energy processes induced by 6 — t- s transitions in the non- 
abelian flavour model predicting pure, CKM-like, left-handed flavour mixing angles in the soft sec- 
tor [32]. In all the plots, all the black points satisfy the constraints of tab. [6] while the red ones 
additionally satisfy BR(i?s — )• < 6 x 10~^. The green points in the plots of BR(ii?s — )• 

vs. BR(i?rf — >• iJL^ iJL~) show the correlation of these observables in the MFV MSSM. 
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lead to the solution of the {g — 2)^ anomaly and large values for the direct CP asymmetry 
^cp(^ si) are typically associated with large values for The departures of S^^ from 

the SM expectations (5*^0)31^1 ~ 0.036 are quite small, in the range < S^^p < 0.1, as expected 
from our model-independent analysis. The corresponding attained values for the semileptonic 
asymmetry normalized to the SM value lie in the range —6 < ^sl/(^sl)sm ^ 6, hence, 
NP effects could be still experimentally visible in A^i^. Yet the confirmation of the large values 
of S^ip observed by CDF and DO will be a serious problem for this kind of models with purely 



left-handed currents. From fig. 14 we also show that models with purely left-handed currents 
can predict very large deviations (in both directions) from the MFV SUSY expectations for 
BR{Bd — )■ yU+/i")/BR(_Bs — > fi'^fi~). In particular, in contrast to the models considered 
so far, the ratio in question can also be significantly larger than in the MFV models and 
BR{Bd — 7- fi~^fi~) can reach values as high as 2 x 10~^, while staying consistent with the 
bound on BR{Bs — )• fi^fi^). 

Finally, in the last row of fig. [T4| we show the correlations involving the CP asymmetries 
in S — )■ iir*/x+/i~, (Ay) vs. (As) and {Aj) vs Sci,Ks- We find a very stringent correlation 
between (A^) and (Ag) with the NP effect in (Aj) being a factor of two larger and with 
opposite sign with respect to (As). This strong correlation is due to the fact that, in the 
considered framework, effects in both (Aj) and (A^) are almost exclusively induced by the 
NP contributions to C7. Moreover, we also find a clean correlation between (Aj) and S(f,Ks 
with —0.2 < {Ai) < —0.1 in the region accounting for the S^Ks anomaly. 

Concerning the hadronic EDMs, the main difference between models with purely LH 
currents with respect to models containing RH currents is that in the former case the quark 
(C)EDMs turn out to be proportional to the external light quark masses. As a result, the 
down- and up-quark (C)EDMs generated by flavour effects are suppressed by small flavour 
mixing angles while not being enhanced by the heaviest quark masses, hence, the hadronic 
EDMs are safely under control. Assuming CKM-like Mis, m = 500 GeV and tan /3 = 10, dd/e 
and are of order ~ 10~^^ cm. The strange quark (C)EDM might reach interesting values 
as it is enhanced, compared to d^/e and d^, by a factor of (ms/m^i) x A^^ ~ 400. However, 
it is still not possible to draw any clear-cut conclusion due to the uncertainties relating dg to 
physical hadronic EDMs. 

In summary, the most striking predictions of the models with purely LH currents as 
opposed to models with large RH currents are 

• the ability to explain the S^Ks anomaly with simultaneous possible explanation of the 
{g — 2)fj_ anomaly and significantly enhanced direct CP asymmetry in 6 — sj, 

• SM-like S^fj,, 

• in the considered SLL model, the possibility of very large deviations for BR(i?rf — >■ 
^'^^~)/BIi{Bs — /U^/U^) in both directions compared to its MFV value, with BR(i?rf — )■ 
fi'^fj'') and BR(i?s — t- fj^^fi^) reaching values as high as 2 x 10~^ and 6 x 10~^, respec- 
tively, 

• Small NP effects in AF = 2 processes, like ex, S^Ks AM^/AM^. Therefore difficulty 
in addressing the UT tension. 
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6.5 Step 5: Comparison with the FBMSSM and the MFV MSSM 



Finally, we want to recall the results for the FBMSSM to discuss the general MFV 
framework and to make a comparison with the models with left-handed CKM-like currents 
discussed in step 4. 

In [68], it was found that the best probes of the FBMSSM are the EDMs of the electron (de) 
and the neutron (dn), as well as flavour changing and CP violating processes in B systems, 
like the CP asymmetries in 6 — )■ 57 and B — )■ (j){ri')Ks, i.e. Acp{b — )■ 57) and S^(^r)')Ks-i 
respectively. The non-standard values for Sif,(^rj')Ksy measured at the B factories, can find a 
natural explanation within the FBMSSM with the effect being typically by a factor of 1.5 
larger in S^Ks m agreement with the pattern observed in the data. 

Interestingly, it was found that the desire of reproducing the observed low values of S(f,Ks 
and Srj'Ks implies: 

• a lower bounds on the electron and neutron EDMs d^^n ^ 10~^^ e cm, 

• a positive and sizable (non-standard) Acp{b — ?■ 57) asymmetry in the ballpark of 1% — 

5%, 

• small NP effects in S^Ks ^^nd AM^/ AMs so that in the FBMSSM these observables can 
be used to extract (3 and \Vtd\, 

• I ex I is enhanced over its SM value at most at a level of < 15%, 

• very small effects in 5*^0 which could, however, be still visible through the semileptonic 
asymmetry A^^, 

• a natural explanation of the Aa^ anomaly under the mild assumption that the slepton 
masses are not much heavier than the squark masses (this is true in almost all the known 
models of SUSY breaking). 

The question we want to address now is, to which extent this pattern of effects gets 
modified within the framework of the MFV MSSM. In particular we will analyse if the richer 
flavour and CP violating structure of the MFV MSSM allows to generate sizeable effects in 



As discussed in sec. 3.1, within a MFV MSSM scenario, the charginos and the charged 
Higgses can induce the FCNC amplitude M^f ^ by means of Ci and C3, with the chargino 



box contributions to C3 given in (3.7) being the only non-negligible contributions sensitive to 
the new phases. 

Moreover, as mentioned in sec. |5} the additional terms in the LL soft masses in the MFV 
MSSM lead to potentially complex Mis 

^ ih + hvDVu , (6.4) 
that will lead to gluino contributions to AF = 2 processes by means of the Wilson coefficient 



Ci given in (3.9) 



In fig. 15, we show the predictions for S^(j, vs. tan /3 in a MFV MSSM scenario taking 
into account the SUSY contributions listed above. As we can see, potentially sizable effects 
to S^^ would be possible only in the very large tan/3 regime (light blue points); however, in 
this case the constraints from both BR.{Bs — )■ /i'''/i~) and BR(6 — )■ 57) become very powerful 
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Figure 15: S"^,^ vs. tan/3 in tlie MFV scenario. Light blue points do not satisfy the S-physics 
constraints while the black ones do. 

and they prevent any visible effect in 5*^^ (black points). This confirms and strengthens the 
general finding of [S^ that within SUSY MFV scenarios, large NP contributions can only be 
expected in AB = 1 transitions. 

The models with left-handed CKM-like currents discussed in step 4 share many similarities 
with the FBMSSM and the MFV MSSM. In particular, all the correlations among 5-physics 
asymmetries are almost the same in these scenarios and also the size of the effects are com- 
parable. As argued in [BH], the {g — 2)^ anomaly is naturally solved once we account for the 
non-standard value for S(f,Ks] this claim has been fully confirmed in the present work where 
we deal with a concrete SUSY breaking scenario, that is a SUGRA scenario. 

The major difference discriminating these scenarios regards their predictions for the lep- 
tonic and hadronic EDMs. Within the FBMSSM, we predicted de,n ^ 10~^^ e cm while solving 
the Sfj)Ks anomaly. On the other hand, within the general class of models with pure left-handed 
CKM-like currents, such a lower bound is significantly relaxed. In fact, in this latter case, 
the source of CP violation is assumed to come from flavour effects and the resulting EDMs 
are suppressed by small flavour mixing angles. This suppression cannot be compensated 
by enhancement factors, as the heavy-light Yukawa ratios, as it happens in presence of RH 
currents. 

Hence, a potential discovery of some of the above CP asymmetries with those peculiar 
correlations without any NP signal in de,n at the level of de,n < 10~^® e cm would most likely 
rule out the FBMSSM and favour non-MFV models with purely left-handed currents. 
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7 Comparison with Other Models 



7.1 Comparison with the LHT and RS Models 

In the present paper we have discussed several specific supersymmetric models that exhibit 
different patterns of flavour and CP violation. In this section we would like to compare briefly 
these patterns with the ones found in the extensive analyses of flavour violation in the Littlest 
Higgs model with T-parity (LHT) [m [721 [731 ESI [761 [77] and in a Randall-Sundrum (RS) 
Model with custodial protection and with a Kaluza-Klein (KK) scale in the reach of the LHC 

[ismn]. 

Let us first recall that the LHT and RS models contain new sources of flavour and CP 
violation but while in the LHT model the operator structure in the effective weak Hamiltonians 
is the same as in the SM [T71 [721 iZS] , in the RS model new operators are present [151 HB] ■ In 
particular the left-right operators contributing to ex imply some fine tuning of the parameters 
in order to keep ex under control if one wants to have the KK scale in the reach of the LHC 
[15] . Analogous tunings, albeit not as serious, are necessary in the RS model, in order to be in 
agreement with the experimental bounds on BR(5 Xs^) [225], BR(/i 67) [226 1 [227 1 [228] 
and EDMs [2291 1230j . that are all dominated by dipole operators. Also contributions to e'/e 
can be large [231j . 

On the other hand in the LHT model BR(i? — )■ Xg'-y) and EDMs are fully under control 
[T7] as NP contributions to the Wilson coefficients of dipole operators are small. Still the 
branching ratio for the decay — )■ 67 can reach the present upper bound [731 EH] and sizable 
contributions to e'/e are possible [3] 

We now compare the predictions of the LHT and RS models with the ones of the super- 
symmetric models analyzed in the present paper. To this end we confine our discussion to 
four very interesting channels: S^^, Bg — )■ — t- tt+z/z/ and yU — j- 67 that are best 
suited for the distinction of supersymmetric models from LHT and RS models. 

• First S^^ can be large in LHT, RS, AC, RVV2 and AKM models, however with the 
following hierarchy for the maximal possible values 

(n \max ^ (Q \raax ^ (O Nmax ^ (O \max ^ / c Nmax (7 ^\ 

Wi/'</>JlHT ~ Wip<l>) AKM ^ WV'0JRVV2 ^ Wip<f>)KS ~ Wip<l>)AC \'-^) 

with (^V'</')lht ~ 0.3 and {S^^)th^ ~ 1 Q 

Instead the values of S^^ in the 5LL model, in the FBMSSM and the MFV MSSM are 
SM-like. 

• While BR(i?s — )■ yU^/i^) in all supersymmetric models discussed by us can be as large 
as the present experimental upper bound, the enhancements of BR(i?s — fj^'^fi') in the 
RS and LHT models do not exceed 10% and 30%, respectively. 

• The opposite pattern is found for — t- tt^uu decays. In all the supersymmetric flavour 
models discussed by us, BIi{K^ — > n^uu) is basically SM-like. On the other hand in 
the RS and LHT models BR(ir+ — )■ n^uu) can be enhanced by as much as a factor 1.5 
and 2.5, respectively. 

^•^In the LHT model, can be significant in spite of tfie purely left-handed structure of the flavour 
violating currents because the 6 — 57 and Bg — !■ ^JtT constraints are easier satisfied than in the (5LL model. 
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• Of particular interest will be the impact of a measurement of S^^ much larger than the 
SM value: 

1. On one hand this would exclude the 5LL, FBMSSM and MFV MSSM models. 

2. On the other hand in the AC and AKM models it would imply a lower bound 
on BR(i?s — )■ fi'^fi~) significantly higher than possible values in the LHT and RS 
models and consequently the measurement of BR(_Bs — could distinguish 
these two classes of models. In the RVV2 model large values of the branching ratio 
in question are possible but not necessarily implied by the anomalously large value 
of S^^. 

3. In the case of — )■ vr+z/z/, a measurement of a sizable S^^ precludes significant 
enhancements of BR(i^'+ — )■ 7r~^i>u) in RS and LHT models. Consequently in such a 
scenario it will be difficult to distinguish these two models from the supersymmetric 
ones, just taking into account the decay — )■ tt^uu. On the other hand, the 
measurement of BR(i^""'" — i- tt+z/z/) much larger than the SM value would favour 
RS and LHT models to all the supersymmetric flavour models discussed by us if 
S^tf) is SM-like but would rule out also the RS and LHT models if S^^ is large. 

• Finally let us recall that while the supersymmetric models and LHT models can reach 
the experimental bound for BR(/i — )■ 67) [73], [76l [77], the ratios BR(yU — )■ 3e)/BR(/i — )■ 67) 
and BR(r — )■ 3/i)/BR(r — )■ 67) are in supersymmetric models roughly by one order of 
magnitude smaller than in the LHT model. 

In summary, we observe that already the combination of S^^f,, BR(i?s — )■ fj^~^fi~) and 
BR(K+ TT+z/z/) should be able to distinguish between the AC, RVV2, AKM, LHT and 
RS models unless all these three observables are found SM-like. The distinction between 
(AC,RVV2,AKM) and (LHT,RS) models can easily be made with the help oiBR{B, 
and BIi{K~^ — t- ti^uu) only but the inclusion in this test of S^^ will be very helpful. The 
distinction between the AC and RVV2 models has been discussed in the previous sections 
while the one between LHT and RS in [121 [77]. Here the correlation between — )■ 7r°z/P 



and — )■ n'^uu markedly different in both models and the hierarchy in (7.1) could play 
important roles. 



7.2 DNA-Flavour Test of New Physics Models 

We have seen in the previous sections and in sec. |7.1| that the patterns of flavour violation 
found in various extensions of the SM differed from model to model, thereby allowing in 
the future to find out which of the models considered by us, if any, can survive the future 
measurements. Undoubtedly, the correlations between various observables that are often 
characteristic for a given model will be of the utmost importance in these tests. 

In tab. [8} we show a summary of the potential size of deviations from the SM results 
allowed for a large number of observables considered in the text, when all existing constraints 
from other observables not listed there are taken into account. We distinguish among: 

• large effects (three red stars), 

• moderate but still visible effects (two blue stars). 
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Table 8: "DNA" of flavour physics effects for the most interesting observables in a selection of SUSY 
and non-SUSY models signals large effects, i^if visible but small effects and implies that 

the given model does not predict sizable effects in that observable. 



• vanishingly small effects (one black star). 

This table can be considered as the collection of the DNA's for various models. These DNA's 
will be modified as new experimental data will be availabe and in certain cases we will be 
able to declare certain models to be disfavoured or even ruled out. 

In constructing the table we did not take into account possible correlations among the 
observables listed there. We have seen that in some models, it is not possible to obtain 
large effects simultaneously for certain pairs or sets of observables and consequently future 
measurements of a few observables considered in tab. |8]will have an impact on the patterns 
shown in this DNA table. It will be interesting to monitor the changes in this table when the 
future experiments will provide new results. 
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8 Flavour vs. LHC Data 



In sec. |6} we have performed a correlated analysis of low-energy predictions, arising in SUSY 
flavour models, for a complete set of flavour observables with the aim of outlining clear 
patterns of deviation from the SM predictions. 

In the following, we will try to address the question of the synergy and interplay existing 
between direct and indirect NP searches. 

For instance, one relevant question could be: which is the expected SUSY spectrum at 
the LHC in case sizable indirect NP effects will be detected in some flavour channels at the 
upcoming experimental facilities? Vice versa, if the LHC will measure the masses of some 
SUSY particles (within a given accuracy), which additional information could we obtain from 
some non-standard effects in flavour observables? 

To answer the above questions, we analyze the flavour physics predictions of the flavour 
models discussed in sec. [6] with respect to the expected spectrum at the LHC. 

In fig. 16, we show the planes for the lightest stop mass (m^J vs. the lightest chargino 
mass (^^xi) ^^^^ — tan /3 planes. The first two rows of fig. 16 correspond to the 

AC model [36] (plots on the left), the RVV2 model [51] (plots in the middle) and the AKM 
model [22] (plots on the right), respectively. The different colors show the possible values 
for 5*^0 in these models as indicated in the overall bar. The last lower plots correspond to 
models with pure left-handed currents with CKM-like mixing angles [12] and the different 
colors indicate the attained values for S^Ks- 

An interesting feature emerging from fig. [16] is that, in the case of the AC, RVV2 and 
AKM models, large effects in S^^ are possible - and indeed even favoured - for a quite 
heavy soft SUSY spectrum, even beyond the LHC reach. However, this is not an evasion of 
the decoupling theorem but, instead, a confirmation of the dominance of the Higgs mediated 
effects in 5^0. In fact, since the effective FCNC couplings of the Higges with the SM fermions 
arise from dimension-4 operators, they do not decouple with the soft SUSY scale. Clearly, the 
double Higgs penguin contributions generating S^^ decouple with the mass of the heavy Higgs 
sector, as confirmed by the second row of fig. 16 , in agreement with the decoupling theorem. 
Moreover, S^^ attains the maximum values for heavy soft masses since the most stringent 
indirect constraints arising from BR(b — )• S7) and ex, decoupling with the soft SUSY sector, 
are relaxed in this case. 

Hence, there exist regions of the SUSY parameter space at the border or even beyond 
the LHC reach where we can expect clear non-standard signals in flavour processes. In these 
regions, flavour phenomena, and thus the indirect search, represent the most powerful tool to 
shed light on NP in case of SUSY. In such a case, where we will obtain only an improved upper 
bound on the scale of SUSY masses, the study of the correlations among flavour observables, 
departing from their SM expectations, would be the only tool at our disposal to unveil the 
nature of the NP theory that is operating. 

An opposite behaviour, compared to the AC, RVV2 and AKM models, is shown by the 
6LL model (last row of fig. 16) where the double Higgs penguins do not play a significant role 
and the flavour observables decouple with the soft SUSY masses. 

Large - non-standard - effects in flavour observables necessarily imply, in this case, a 
SUSY spectrum within the LHC reach. Hence, combining the information on the spectrum 
from the LHC with the information from the low-energy flavour processes, we would be 
able, to some extent, to measure the mixing angle regulating the flavour transition. Such 
an achievement would represent a crucial step forward towards the reconstruction of the 
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underlying flavour symmetry at work. The described situation represents one clear example 
of the synergy and interplay of direct and indirect NP searches. 



9 Summary 

The coming years will witness tremendous progress at the high energy frontier accomplished 
primarily at the LHC, where the available energy will be increased by one order of magnitude. 
Equivalently, for the first time we will be able to resolve directly distances well below 10~^^ m, 
that have been explored so far. Parallel to these developments, important advances are 
expected at the high precision frontier through the improved i?s-physics experiments at the 
Tevatron and in particular LHCb at CERN. At later stages in the coming decade these efforts 
will be strengthened by new rare K experiments at J-PARC, the NA62 collaboration at CERN 
and possibly Project X at Fermilab as well as Belle-II at KEK and the planned Superi? facility 
in Rome. The latter two will also provide new insights into the FCNC processes in the D 
meson system and in charged lepton decays. 

While the main goal at the high energy frontier is the discovery of new particles and the 
determination of their masses, the main goal of flavour physics is the search for the footprints 
of these new particles in rare processes and the determination of their couplings. As the latter 
exploration of very short distance scales is indirect, only measurements of a large number of 
observables and the study of correlations between them in a given extension of the SM and 
in particular of patterns of flavour violation characteristic for a given model can allow us to 
identify the correct NP scenario. 

In the present paper, we have performed an extensive study of processes governed by 6 — s 
transitions and of their correlations with processes governed by 6 — ?• c? transitions, s ^ d 
transitions, — oscillations, lepton flavour violating decays, electric dipole moments and 
{g — 2)^. To this end, we have considered a number of supersymmetric flavour models that 
on one hand aim at the explanation of the observed hierarchical fermion masses and mixings 
and on the other hand provide natural suppression of FCNC transitions. In particular we 
have analyzed the following representative scenarios: 

i) dominance of right-handed currents (abelian model by Agashe and Carone [3^), 

ii) comparable left- and right-handed currents with CKM-like mixing angles represented 
by the special version (RVV2) of the non abelian SU (3) model by Ross, Velasco and 
Vives |51j as discussed recently in [62j and the model by Antusch, King and Malinsky 
(AKM) 



iii) dominance of left-handed currents in non-abelian models 

In the choice of these three classes of flavour models, we were guided by our model indepen- 
dent analysis, as these three scenarios predicting quite different patterns of flavour violation 
should give a good representation of most SUSY flavour models discussed in the hterature. 

The distinct patterns of flavour violation found in each scenario have been listed in detail 



in sec. 6 and the corresponding plots can be found in flgures 11-14 



The main messages from our analysis of the models in question are as follows: 

• Supersymmetric models with right-handed currents (AC, RVV2, AKM) and those with 
left-handed currents can be globally distinguished by the values of the CP-asymmetries 
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Figure 16: Tlie plane of the lightest stop mass {mf_^) vs. the lightest chargino mass (m^j^j) as well 



as the Mh — tan /3 planes. The first two rows of fig. 16 correspond to the AC model [36] (plots on 
the left), the RVV2 model [M] (plots in the middle) and the AKM model [52] (plots on the right), 
respectively. The different colors show the possible values for S*^^ in these models as indicated in the 
overall bar. The last lower plots correspond to models with pure left-handed currents with CKM-like 
mixing angles [12] and the different colors indicate the attained values for S(j,Ks- 



68 



S^(j) and S^Ks with the following important result: none of the models considered by 
us can simultaneously solve the S^^ and S^Ks anomalies observed in the data. In the 
models with RH currents, S^^ can naturally be much larger than its SM value, while 
S^Ks remains either SM-like or its correlation with S^^ is inconsistent with the data. 
On the contrary, in the models with LH currents, S*^^ remains SM-like, while the S'^i^-g 
anomaly can easily be solved. Thus, already future precise measurements of S^^ and 
S^Ks will select one of these two classes of models, if any. 

The desire to explain the anomaly within the models with RH currents unambigu- 
ously implies, in the case of the AC and the AKM models, values of BR(i?s as 
high as several 10~^, while in the RW2 model, such values are also possible but not nec- 
essarily implied by the large value of S'^^. Interestingly enough, in all these models large 
values of S^^ can also explain the {g — 2)^ anomaly. Moreover, in the AC and RVV2 
models, the ratios BR(_Brf -)> )/BR(_Bs /U^/i") and BR{Bd fi'^fi^)/AMs 

can significantly differ from the MFV prediction, providing a splendid opportunity 
to shed light on new sources of flavour violation beyond the CKM ones. Further, 
BR{Ba — > Ijl'^Ijl~)/BR{Bs — > jJ'^IJr) is predicted to be dominantly smaller than its 
MFV value especially for BR(5, /x+zx") > 10"^. Values of BR{Ba n^n') as high 
as 1 X 10~^ are still possible in these models. 

The hadronic EDMs represent very sensitive probes of SUSY flavour models with right- 
handed currents. In the AC model, large values for the neutron EDM might be easily 
generated by both the up- and strange-quark (C)EDM. In the former case, visible CPV 
effects in — mixing are also expected while in the latter case large CPV effects 
in Bg systems are unavoidable. The RVV2 and AKM models predict values for the 
down-quark (C)EDM and, hence for the neutron EDM, above the ^ 10~^^ecm level 
when a large - non-standard - 5"^^ is generated. All the above models predict a large 
strange-quark (C)EDM, hence, a reliable knowledge of its contribution to the hadronic 
EDMs, by means of lattice QCD techniques, would be of the utmost importance to 
probe or to falsify flavour models embedded in a SUSY framework. 

In the RVV2 and AKM models, the desire to explain the {g — 2)^ anomaly implies 
that BR(/x — )■ 67) > 10~^^, in the reach of the MEG experiment, and de > 10~^^ ecm 
{de < 10~^°ecm) in the RW2 (AKM) model is predicted. Moreover, in the case of 
the RVV2 model, BR(t /X7) > 10~^ is then in the reach of Super-B machines, while 
T ^ fi'-f remains most likely beyond the Super-B reach in the case of the AKM model. 
The explanation of the {g — 2)^ anomaly, combined with non-standard values for S^,f,, 
would imply BR(;U 67) > 10~^^ in both the RVV2 and AKM models. 

Finally, while the abelian AC model resolves the present UT tension through the modifi- 
cation of the ratio AMd/AM,^, the non-abelian flavour models RVV2 and AKM provide 
the solution through NP contributions to ex- Moreover, while the AC model predicts 
sizable NP contributions to D'^ — mixing, such contributions are negligibly small in 
the RW2 and AKM models. 

In the supersymmetric models with LH currents only, the desire to explain the S^f^Ks 
anomaly unambiguously implies that the direct CP asymmetry in 6 — > 57 is much 
larger than its SM value. This is in contrast to the models with RH currents where this 
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asymmetry remains SM-like. Also the solution to the {g — 2)^ anomaly can be easily 
accounted for. 

• Interestingly, also in the LH-current-models -that are in many aspects similar to MFV 
models-, the ratio BR(i?rf — )■ /i"'"/i~)/BR(i?s — )■ can not only deviate significantly 
from its MFV value, but in contrast to the models with RH currents considered by us 
can also be much larger than the latter value. Consequently, values for 'BK{B(i — )■ fi~^fi~) 
as high as 2 X 10~^ are still possible while being consistent with the bounds on all other 
observables, in particular the one on BR(i?s — ^^^~). Also interesting correlations 
between S^Ks CP asymmetries in B ^ K*i'^£~ are found. 

• The FBMSSM as well as the MFV MSSM bear several similarities to the models with 
LH currents only. In particular, in all these models, S^rp is basically SM like. The 
FBMSSM and the MFV MSSM differ however from the 6LL model especially in their 
predictions of the EDMs. In the FBMSSM and the MFV MSSM, large effects in CP 
violating AF = 1 observables naturally imply values for the EDMs close to the current 
experimental limits. In the 6LL model, on the other hand, the EDMs can be easily well 
under control. 

This summary of our most important results already shows that the simultaneous study 
of various flavour violating processes can allow us to distinguish various NP scenarios. More 
results can be found in the numerous flgures presented by us. 

We have also presented a comprehensive model-independent analysis of the MSSM with 
general sources of flavour and CP violation that goes in several aspects, listed at the beginning 
of our paper, beyond many analyses found in the literature. The bounds on the mass insertions 



i^u,d)i,j be found in flgures [g]- 9 



Any version of the MSSM based on a SUGRA spectrum 
has to satisfy them. Beyond our MSSM analysis, we have also emphasized the usefulness of 
the -Rfc — 7 plane in exhibiting transparently various tensions in present UT analyses. 

It will be exciting to monitor upcoming results from Tevatron and the LHC on S^^ and 
Bg — . Already these two measurements will be capable of excluding some of the models 
considered by us and distinguish them from the LHT and RS model with custodial protection 
in which 5*^^ can be large but Bg^d — ^ Z^"*"/^" remain SM-like. Further observables analyzed 
by us will help to identify more precisely the correct extension of the SM. In particular, while 
the branching ratios for K — > ttui? decays in the supersymmetric models considered by us 
remain SM-like, they can be signiflcantly enhanced in the LHT and RS models. 

A DNA-Flavour Test proposed by us should give still a deeper insight into the patterns of 
flavour violation in various scenarios, in particular when it is considered simultaneously with 
various correlations present in concrete models. The interplay of these efforts with the direct 
searches for NP will be most exciting. 

In conclusion, the above physics cases are representative of the richness which is present 
in flavour physics once we assume flavour models embedded within a gravity mediated SUSY 
breaking framework. It could be that, at the end, flavour physics is one of the very few tools 
we have to understand from low-energy physics whether nature possesses supersymmetry at 
the root of its symmetries, or whether other scenarios like LHT or RS are realized in nature. 
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A Explicit Expressions for the Loop Functions 

Loop Functions for the AF = 2 Mixing Amplitudes 

(1) , , 11 + 144x + 27x2 - 2x3 3 ^ 

q) X) = z -; z r^logX , A.l 

^ ' 108(1 -x)4 18(1 -x)5 ^ ' ^ 

(2) , , 33 + 665x + 237x2 - 39x3 ^ 4^4 ^(^26 + 49x) , 

fifi (x) = z h logx , (A. 2 

^ ' 216(1 -x)5 18(1 -x)6 ^ ' ^ 

(3) , , 66 + 1835x + 1005x2 - 255x=^ + 55x^ - 6x^ x(13 + 32x) , , , ^ 

g\-'>(x) = ^— ^logx, (A. 3 

^ ' 1080(1 -x)6 18(1 -x)7 ^ ' ^ 

(i), , 2 - 99x - 54x2 + 7x3 _^ iq^\ 

Qd (x) = } rz r/logX, (A.4 

^ ^ 18(1 -x)4 3(1 -x)5 ^ ' ^ 

(2) , , 3 - 212x - 192x2 + 48x3 - 7x^ 10x(l 5x) , 

= + 3(1 -x)6 ' ' ^^-^ 

(3) , , 12 - 1117x - 1452x2 + 528x3 - 152x^ + 21x^ x(5 + 31x) , 

9a (x) = z \ r/logx, (A. 6 

^ ' 180(1 -x)6 3(1 -x)7 ^ ' ^ 

(1),, 10 + 117x + 18x2 - x3 x(ll + 13x), 

= W^' -w^'°<'^- <*-^ 

(2),, 15 + 272x + 84x2 - 12x3 + x^ 2x(ll + 19x) , 

9^ (x) = z 1 z — logx , A. 8 

^ ^ 54(1 -x)5 9(1 -x)6 ^ ' ^ 

(3) , , 60 + 1507x + 732x2 - 168x3 + 32x^ - 3x^ x(ll + 25x) , , , ^ 

"^f^)" sw^^ 9(1 

, , , 4(l + x) 8x , 

hJx) = ^ i7z + — rrrlogX, (A.IO 

^ ' 3(1 -x)2 3(1 -x)3 ^ ' ^ 
, , , 4(2 + 5x-x2) 8x , 

= "2(1^- 2(13^ '"i^"' 
^> = - (1 - x)(l - ,) + (1 - xA, - X) ^ + (1 - - ,) • <*-^3 
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, . X+1 X 

^ , , x'^ -6x-17 3x + 1 , 
^^^"^^ 6(1-.)^ -(1^^°^-' 



Loop Functions for b ^ sj 



^^^"^ = -12(13^- 6(1^ 

= -4(13^ + 2(13^^°^^' 
^(2), , 13 -7a; 3 + 2x - 2a;2 

^"^ = -24(1^- 12(1-.)^ 

(2) , , 1 + 5a; a;(2 + x) 

/S(^,l/) = -^(/?,8(^)-/?8(l/)) , 

(1), , 2(1 + 5x) 4a;(2 + x) , 
("^ = -9(1^-9(1-^^°^"' 

(2), , 2(1 + 10x + ,t2) 4a;(l + a;), 

= 9(1-.)^ - 3(1^ ' 

m, , 11 + ,T 9 + 16,T-a;\ 
^« ("^ = 3(1^+ 6(1 -x)^ 
(2), , 53 + 44a;-a;2 3 + Ha; + 2a;2 , 

= 12(1 -x)^ + 2(1 -x)^ ■ 

Loop Functions for ii ij^ 

/2„(a;) = (-5a;2 + 4a; + 1 + 2a;(a; + 2) loga;)/(4(l - xf) , 
/3„(a;) = (1 + 9x - 9a;^ - + 6a;(x + 1) loga;)/(3(l - xf) , 
f2c{x) = (-x' - 4a; + 5 + 2(2x + 1) logx)/(2(l - x)^) . 
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Loop Functions for the EDMs 



S "(3)^^N 1 f(3)/^N ^ 3 (3) 



6^ 



(o: 



.(1 

Jo 



(2: 



f(3 

Jo 



(o: 



n 



X 



X 



(x 



X 



X 



X 



1 — x^ + 2x log X 
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-1 - 4x + - 2x(x + 2) log X 
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-10 - 9,r + 18.r- + .r'^ - 3(1 + 6,r + 3,r-) logx 
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SM Loop Functions 
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3x^ log X 
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+ 
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logx 
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B MSSM parameter convention 



Throughout sec. [3] we adopted the convention of |ll3j for squark matrices and trihnear terms. 
In tab. 9 we show a dictionary between this convention and the "SUSY Les Houches Accord 
2" (SLHA 2) convention in [2^2] . 



SLHA 2 [232J 


m 


Tu, Td, Te 

m^, TTL j 

ml,m\, ml 


~^Z-> +^1; 

{mlf,{mlf,{mlY 

{Mir, {Mir 



Table 9: Dictionary between the SLHA 2 convention of |232j and the convention of |113] 
Note that At and Ab used in sec. j3jare defined as mtAt = ^ {^u)s3 and mbAb = ^ (^d)33- 
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